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Abstract 

Statistical systems displaying a strongly anisotropic or dynamical scaling behaviour are char- 
acterized by an anisotropy exponent 9 or a dynamical exponent z. For a given value of 9 (or 
z), we construct local scale transformations, which can be viewed as scale transformations with 
a space-time-dependent dilatation factor. Two distinct types of local scale transformations are 
found. The first type may describe strongly anisotropic scaling of static systems with a given 
value of 9, whereas the second type may describe dynamical scaling with a dynamical exponent 
z. Local scale transformations act as a dynamical symmetry group of certain non-local free-field 
theories. Known special cases of local scale invariance are conformal invariance for 9 = 1 and 
Schrodinger invariance for 9 = 2. 

The hypothesis of local scale invariance implies that two-point functions of quasiprimary 
operators satisfy certain linear fractional differential equations, which are constructed from 
commuting fractional derivatives. The explicit solution of these yields exact expressions for 
two-point correlators at equilibrium and for two-point response functions out of equilibrium. A 
particularly simple and general form is found for the two-time autoresponse function. These 
predictions are explicitly confirmed at the uniaxial Lifshitz points in the ANNNI and ANNNS 
models and in the aging behaviour of simple ferromagnets such as the kinetic Glauber-Ising 
model and the kinetic spherical model with a non-conserved order parameter undergoing either 
phase-ordering kinetics or non-equilibrium critical dynamics. 
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1 Introduction 



Critical phenomena are known at least since 1822 when Cagnard de la Tour observed critical opales- 
cence in binary mixtures of alcohol and water. The current understanding of (isotropic equilibrium) 
critical phenomena, see e.g. |37 , |116| , |34| , H| , is based on the covariance of the n-point correlators 



G n = G(ri, . . . , r n ) = (0 1 (r 1 ) . . . (f) n (r n )) under global scaling transformations r< — > bVi 

G(b ri , br n ) = b-^+-+ x ")G( ri , . . . , r n ) (1.1) 

precisely at the critical point. Here the fa are scaling operators^ with scaling dimension Xi and one 
might consider the formal covariance of the fa 

fa{br t ) = b-^iri) (1.2) 

as a compact way to express the covariance fll.ip of the correlators. In isotropic (e.g. rotation- 
invariant) equilibrium systems, the fa corresponds to the physical order parameter or energy den- 
sities and so on. Eq. (|1.1|) may be derived from the renormalization group and in turn implies the 
phenomeno logical scaling behaviour of the various observables of interest. It has been known since 
a long time that in systems with sufficiently short-ranged interactions, G n actually transform co- 
variantly under the conformal group, that is under space-dependent or local scale transformations 
b = b(r) such that the angles are kept unchanged |92]| . Since in two dimensions, the Lie algebra of 
the conformal group is the infinite-dimensional Virasoro algebra, strong constraints on the possible 
fa present in a 2D conformally invariant theory follow ||. Roughly, for a given unitary 2D conformal 
theory, the entire set of scaling operators fa, the values of their scaling dimensions Xi and the critical 
n-point correlators G n can be found exactly. Furthermore, there is a classification of the modular 
invariant partition functions of unitary models at criticality (the ADE classification) which goes a 
long way towards a classification of the universality classes of 2D conformally invariant critical points 
(for reviews, see |0J ^T], |57|). 



Here we are interested in critical systems where the n-point functions satisfy an anisotropic scale 
covariance of the form 

G(b%, br <!,..., b%, br n ) = b-^ + -^G{t x , r u ...,t n , r n ) (1.3) 

where we distinguish so-called 'spatial' coordinates r*j and 'temporal' coordinates ij. The exponent 



6 is the anisotropy exponent. By definition, a system whose n-point functions satisfy (1.3) with 
9 7^ 1 is a strongly anisotropic critical system. Systems of this kind are quite common. For example, 
eq. ( |1.3|) is realized in (i) static equilibrium critical behaviour in anisotropic systems such as dipolar- 
coupled uniaxial ferromagnets |2| and/or at a Lifshitz point in systems with competing interactions 
|64|, p.07| , pi)] or even anisotropic surface-induced disorder ||110|| , (ii) anisotropic criticality in steady 



states of non-equilibrium systems such as driven diffusive systems ||103 , 77], stochastic surface growth 
72| or such as directed percolation. In these cases, r and t are merely labels for different directions 
in space and the G n are in the case (i) equal to the n-point correlators C n of the physical scaling 
operators, (iii) Further examples are found in quantum critical points, see |98|. Anisotropic 
scaling also occurs in (iv) critical dynamics of statistical systems at equilibrium [[51] or (v) in non- 
equilibrium dynamical scaling phenomena [JIBL [93], [77], |102| , |2Uj. In the cases (iv) and (v), t 



2 We follow the terminology of ]l9| ]: the (f>i are called scaling operators, because if the theory is quantized in the 
operator formalism, (pi — > (f>i becomes a field operator. The variables hi canonically conjugate to (pi are called the 
scaling fields. For the Ising model, the scaling operator (p a = a is the order parameter density and its canonically 
conjugate scaling field h is the magnetic field. 
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represents the physical time and the system's behaviour is characterized jointly by the time-dependent 
rz-point correlators C n as well as with the n-point response functions R n . Habitually, 9 = z then is 
referred to as the dynamical exponent. At equilibrium, the C n and R n are related by the fluctuation- 
dissipation theorem |7T], [HJ , but no such relation is known to hold for systems far from equilibrium. 

We ask: is it possible to extend the dynamical scaling ( |i.3| ) towards space-time-dependent rescal- 
ing factors b = b(t, r) such that the n-point functions G n still transform covariantly ? 

It is part of the problem to establish what kind(s) of space-time transformations might be sen- 
sibly included into the set of generalized scaling transformations. Also, for non-static systems, 
covariance under a larger scaling group may or may not hold simultaneously for correlators and 
response functions. Another aspect of the problem is best illustrated for the two-point function 
G*2 = G(ti,ri,t2,r2) = G(t,r) where for simplicity we assume for the moment space and time 
translation invariance and t — t\ — ti and r = r% — r%. From (|1.3|) one has the scaling form 

G 2 = G(t, r) = t- 2x ' e g{u) , u = \r\ e /t (1.4) 

where - in contrast to the situation of isotropic equilibrium points with 9 = 1 (see below) - the scaling 
function Q{u) is undetermined. We look for general arguments which would allow us to determine 
the form of Q, independently of any specific model. In turn, once we have found some sufficiently 
general local scaling transformations, and thus predicted the form of G 2 , the explicit comparison 
with model results, either analytical or numerical, will provide important tests. Several examples of 
this kind will be discussed in this paper. 



Some time ago, Cardy had discussed the presence of local scaling for critical dynamics [fTS 
Starting from the observation that static 2D critical systems are conformally invariant, he argued that 
the response functions should transform covariantly under the set of transformations r — > b(r)r and 
t — > b(r) e t. Through a conformal transformation, the response function was mapped from 2D infinite 
space onto the strip geometry and found there through van Hove theory. Explicit expressions for 
the scaling function Q were obtained for both non-conserved (then Q(u) ~ e~ u , up to normalization 



constants) and conserved order parameters ||18|| . However, these forms have to the best of our 
knowledge so far never been reproduced in any model beyond simple mean field (i.e. van Hove) 
theory. That had triggered us to try to study the construction of groups of local anisotropic scale 
transformations somewhat more systematically, beginning with the simplest case of Schrodinger 



invariance which holds for 9 = 2 [52, B3f]. At the time, it appeared to be suggestive that the exactly 



known Green's function of the ID kinetic Ising model with Glauber dynamics |43[] could be recovered 



this way. How these initial results might be extended beyond the 9 = 2 case is the subject of this 
paper. 

The outline of the paper is as follows: in section 2, we shall review some basic results of conformal 
invariance and of the simplest case of strongly anisotropic scaling, which occurs if 9 = 2. In this case, 
there does exist a Lie group of local scale transformations, which is known as the Schrodinger group 
54] , Building on the analogy with this case and conformal invariance for 9 = 1, we discuss in 



section 3 the systematic construction of infinitesimal local scale transformations which are compatible 



with the anisotropic scaling (O). We shall see that there are two distinct solutions, one corresponding 
to strongly anisotropic scaling at equilibrium and the other corresponding to dynamical scaling. 
We also show that the local scale transformations so constructed act as dynamical symmetries on 
some linear field equations of fractional order. Furthermore, linear fractional differential equations 
which are satisfied by the two-point scaling functions Q{u) are derived. In section 4, these are 
solved explicitly and the form of Q{u) is thus determined. In section 5, these explicit expressions are 
tested by comparing them with results from several distinct models with strongly anisotropic scaling, 
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notably Lifshitz points in systems with competing interactions such as the ANNNI model and for 
some ferromagnetic non-equilibrium spin systems (especially the Glauber-Ising model in 2D and 3D) 
undergoing aging after being quenched from some disordered initial state to a temperature at or below 
criticality. We also comment on equilibrium critical dynamics. A reader mainly interested in the 
applications may start reading this section first and refer back to the earlier ones if necessary. Section 
6 gives our conclusions. Several technical points are discussed in the appendices. In appendix A we 
discuss the construction of commuting fractional derivatives and prove several simple rules useful for 
practical calculations. In appendix B we generalize the generators of the Schrodinger Lie algebra to 
d > 1 space dimensions and in appendix C we present an alternative route towards the construction of 
local scale transformations. Appendix D discusses further the solution of fractional-order differential 
equations through series methods. 



2 Conformal and Schrodinger invariance 

Our objective will be the systematic construction of infinitesimal local scale transformations with 
anisotropy exponents 6^1. Consider the scaling of a two-point function 

G = G(t, r) = b 2x G(b e t, br) = r 2x/e $(rr l/e ) = r - 2x {l(tr- e ) (2.1) 

where t = t\ — t 2 , f = fi — ^2, t = \r\ and a; is a scaling dimension. For convenience, we also assumed 
spatio-temporal translation invariance. The scaling of G is described by the scaling functions $(m) or 
alternatively by f2(t>) = v~ 2x / e & (f -1 / 61 ). In this section we concentrate on the formal consequences 
of the scaling Q2.1D and postpone the question of the physical meaning of G to a later stage. 

Considering ( |2.1| ) for r = 0, one has G(t, 0) ~ f~ 2x / e anc i if £ = 0, one has G(0,r) ~ r~ 2x . 
Therefore, 

$(0) = $ , $(m) ~ &ooU~ 2x ; u — > 00 

fi(0) = n , Q(v) ~ ftootT 2x/e ; v -> 00 (2.2) 

where $0,00 and Qq i00 are generically non-vanishing constants. This exhausts the information scale 
invariance alone can provide. 



2.1 Conformal transformations 

Consider static and isotropic systems with short-ranged interactions. Then the two-point function 
G(t,r) = (0i (ti, ri)(j) 2 (t2, r 2 )) is the correlation function of the physical scaling operators 4>i,i- If 
these (pi are actually quasiprimary || scaling operators, G does transform covariantly under the 
action of the conformal group. To be specific, we restrict ourselves to two dimensions (here t and r 
merely label the different directions) and introduce the complex variables 

z = t + ir , z = t — ir (2.3) 

Then the projective conformal transformations are given by 

z ^ z ' = ^±Jl . a s-f3 1= i (2.4) 
7Z + 
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and similarly for z. Writing z' = z'(z) = z + £{z), the infinitesimal generators read 

4 = -z n+l d z , 4 = (2.5) 
and satisfy the commutation relations 

[4j ^m] = ( n ~ m )£n+m j [Ai> ^m] = , = ( n ~~ m )£n+m (2-6) 

In fact, although the generators £ n were initially only constructed for n = —1,0,1, the £ n can be 
written down for ailneZ and ( |2.6|) still holds. The existence of this infinite-dimensional Lie algebra, 
known as the Virasoro algebra without central charge, is peculiar to two spatial dimensions. The set 
( p.4| ) corresponds to the finite-dimensional subalgebra {4ti,o> 4ti,o}- 

The simplest possible way scaling operators can transform under the set ( |2.4|) is realized by the 
quasiprimary operators [[LOO] , || , which transform as 



6</>i(z, z) = (Aie'iz) + e(z)d z + Aie'iz) + e(z)d z ) z) (2.7) 

where A, and Aj are called the conformal weights of the operator fa. If <pi is a scalar under (space- 
time) rotations (we shall always assume this to be the case), Aj = Aj = Xi/2, where Xi is the scaling 
dimension of 4>i- If z(z) = £ z n+1 , one then has 5<fii(z, z) = —e[£ n + £ n )<f>i(z, z) where the generators 
£ n , I n now read 

£ n = -z n+1 d z - A,(n + l)z n , £ n = -z n+l d- z - A,(n + l)z n (2.8) 

and again satisfy (|2.6|) . Later on, we shall work with the generators 

X n . £ n ~\~ £ n , Y n . i(^n £n) (2.9) 

which satisfy the commutation relations 

[X n , X m ) = (n - m)X n+m , [X n , Y m ] = (n - m)Y n+m , [Y n , Y m ] = -(n- m)X n+m (2.10) 

The covariance of G under finite projective conformal transformations leads to the projective con- 
formal Ward identities for the n-point functions G of quasiprimary scaling operators (see [|100|| for a 
detailed discussion on quasiprimary operators) 

£ n G = £ n G = < — ► X n G = Y n G = (2.11) 

for n = ±1, and the generators as defined in eqs. (|2.8|J2.9| ). This gives for the two-point function of 
two scalar quasiprimary operators |92] 

G(M2;ri,r a ) = G 12 5 X1>X2 {{ Zl - z 2 )(z 1 - z 2 )y Xl = G 12 5 XljX2 ((h - t 2 f + (n - r 2 ) 2 )^ (2.12) 

where G\ 2 is a normalization constant (usually, one sets G\ 2 = 1). Comparison with ( p.l|) gives the 
scaling function fl(v) ~ (1 + v 2 )~ x . The constraint x\ = x 2 is the only result which goes beyond 
simple scale and rotation invariance. 

The three-point function is |)2 



G(ti,t 2 ,t 3 ;n,r 2 ,r 3 ) = G 123 P^ 123 P^ 231 P3? 312 (2.13) 

where 

p 2 ab = z ab z ab = (t a - t b ) 2 + (r a - r b f , z ab = z a -z b (2.14) 
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and x a b c := x a + x^ — x c . The constant G123 is the operator product expansion coefficient of the three 
quasiprimary operators 0i,2,3- For scalar quasiprimary operators, the results ( |2.12|j2.1^ ) remain also 
valid in d > 2 dimensions, since two or three points can by translations and/or rotations always be 
brought into any predetermined plane. 

The conformal invariance of scale- and rotation-invariant systems is well established. A convenient 
way to show this proceeds via the derivation of Ward identities, invoking the (improved) energy- 
momentum tensor. These Ward identities hold for systems with local interactions and it can be 
shown that any n-point function which is translation-, rotation- and scale invariant is automatically 
invariant under any projective conformal transformation, see e.g. |19 , [34] , 41 , 57fl . 

We have restricted ourselves to quasiprimary operators || , which is all what we shall need in this 
paper. 



2.2 Schrodinger transformations 



50p by the following set of trans- 
a<y-/?7 = l (2.15) 



The Schrodinger group in d + 1 dimensions is usually defined 
formations 

. TZr + vt + a . at + P 

r — > r = , t—>t= - 

jt + 5 -ft + 5 

where a, (3, 7, 5, v, a are real parameters and TZ is a rotation matrix in d spatial dimensions. The 
Schrodinger group can be obtained as a semi-direct product of the Galilei group with the group 
Sl(2, R) of the real projective transformations in time. A faithful d + 2- dimensional matrix represen- 
tation is 

TZ v a 



a (3 
7 5 



c 



yy 



(2.16) 



According to Niederer [83], the group ( 2.15 ) is the largest group which transforms any solution of 
the free Schrodinger equation 

d 1 d d . 
1— + — — •— 1^ = (2.17) 



dt 2m dr dr t 
into another solution of ( [2.17|) through (t, r) 1— > g(t, r), ip hh> T g ip 

(Tgij) (t,r) = f g (g- 1 (t,r))ij(g-\t,r)) 



where 84, S 7 



f g (t,r) = {jt + S) 



-d/2 



exp 



im 7r 2 + ITZr ■ (7a — 5v) + 7a 2 — t5v 2 + 2-fav 

'1 rt+~s 



(2.18) 



(2.19) 



Independently, it was shown by Hagen [|5(J that the non-relativistic free field theory is Schrodinger- 
invariant (see also |78[). Furthermore, according to Barut the Schrodinger group in d space 
dimension can be obtained by a group contraction (where the speed of light c — > 00) from the 
conformal group in d+1 dimensions (this implies a certain rescaling of the mass as well). Formally, 
one may go over to the diffusion equation by letting m = (2iD)~ 1 , where D is the diffusion constant. 

In order to implement the Galilei invariance of the free Schrodinger equation and of a statistical 
system described by it, the wave function ip and the scaling operators (pi of such a theory will under 
a Galilei transformation pick up a complex phase as described by f g 7^ and characterized by the 
mass m ||, |75| . By analogy with conformal invariance M, we call those scaling operators with the 
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simplest possible transformation behaviour under infinitesimal transformation quasiprimary, that is 
8x4>i = —zXnifii and 6y4>i = —eY m <f>i. In d = 1 space dimensions, to which we restrict here for 



simplicity (then 1Z = 1), we have [[53 
X„ 



Y 

- 1 it. 



-t n+1 d t - 



-t m+1 / 2 d r 

-Mt n 



n + 1 _ _ n(n + 1) . . _ 9 x, _. r 

t n r<9 r i '-Mt^r 2 (n + l)e 

2 4 2 



m + - ] .Mt m ~ 1/2 r 



^2.20) 



for a quasiprimary operator with scaling dimension x and 'mass 7 M. = im. Here x and Ai are 
quantum numbers which can be used to characterize the scaling operator <p. Extensions to spatial 
dimensions d > 1 are briefly described in appendix B. Necessarily, any Schrodinger-invariant theory 
contains along with also the conjugate scaling operator 0*, characterized by the pair (x, —M). For 
x = M. = 0, we recover the infinitesimal transformations of the Lie group ( |2.15|) . The commutation 
relations are 



(n - m)X n+ri 
(n — m)M n+r 



[Y n , M„ 



- - m ) Y n+m , [X n , M„ 



-mM„ 



n+m 



[M n ,M m ] = 



(2.21] 



in d = 1 space dimensions. The infinitesimal generators of the finite transformations ( [2.15 ) are given 
by the set {^±i,o> ^±i/2> M }. In this case the generator M commutes with the entire algebra. The 
eigenvalue — M. of M can be used along with the eigenvalue Q of the quadratic Casimir operator 



87 



Q := (aM X - 2{Y_i , Yi }) - 2 ^2M X^ - Y\,2M X 1 - Y|} (2.22) 
(where {A, B} := AB + BA) to characterize the unitary irreducible (projective) representations of 



the Lie algebra ( p.21|) of the Schrodinger group ||87|| . One can show that the representations with 
Q = realized on scalar functions reproduce the transformation ( 2.18 . 2~T9| ) [^7]. Frequently, the 
algebra with M n ^ is referred to as a centrally extended algebra. However, since the algebra (p.21[) 
with M n = (i.e. M. = 0) is not semi-simple, its central extension is quite different from those of the 
conformal algebra (|2.6|) . Since for the physical applications, we shall need M. ^ anyway, we shall 



refer to (|2.21| ) as the Schrodinger Lie algebra tout court and avoid talking of any 'central extensions' 
in this context. 

As was the case for 2D conformal transformations, one may write down the generators X n , M n 



for any fiGZ and Y m for any m E Z + \ such that ( |2.21| ) remains valid |52|, |53| . 

By definition ||53|| , n-point functions G of quasiprimary scaling operators <pi with respect to the 
Schrodinger group satisfy 

X n G = Y rn G = (2.23) 

with n = —1, 0, 1 and m = —1/2, +1/2. Consequently, the only non-vanishing two-point function of 
scalar quasiprimary scaling operators is, for any spatial dimension d > 1, 



hih, rx)(f)* 2 (t2, r 2 )) = G 12 6 Xl>X2 5. 



M1M2 



{ti-h)- xi exp 



Mx (r-i - r 2 ) 



U > to 



(2.24) 



(4>*4>*) = provided Mi ^ |53l . Usually, the normalization constant Gu = 1. 

~ e~ Mu I 2 . Similarly, the basic 



whereas 

Comparison with the form (|2.1|) gives the scaling function $(n) 
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non-vanishing three-point function of quasiprimary operators reads 



.(h, rx)4> 2 {t 2 , r 2 )(f>l{h, r 3 )> = S Ma+MbMc t 13 J 



X132/2 .—1231/2 .-1123/2 



^23 



12 



exp 



All rf 3 A4 2 T*| 3 



2 t 



13 



2 i 



23 



'12,3 



Tl3*23 _ r 23^ 



13 J 



^12^13^23 



ti > , to > iq 



(2.25) 



+ x b — x c and Q a b,c is an arbitrary differentiable scaling 
'(ft*), while (000) = (0*0*0*) = 0, unless the 'mass' M 



with t a6 = t a - t b) r ab = r a - r b) x abc = x 
function. A similar expression holds for 
of the scaling operator vanishes. 

It is instructive to compare the form of the two- and three-point functions ( |2.12| , pTT3| ) as obtained 
from conformal invariance with the expressions ( p.24| , pT25l) following from Schrodinger invariance. 
As might have been anticipated from comparing the finite transformations fl2.4| ) and ( |2.15| ), the 
dependence on z ab and t ab , respectively, is identical provided the scaling dimensions Xi are replaced 
by Xi/6. For the two-point function, we have in both cases the constraint x a = x b and one could 
extend the arguments of ||100|| on n-point functions between derivatives of quasiprimary operators 
from conformal to Schrodinger invariance. On the other hand, Schrodinger invariance yields the 
constraints M. a = M. b for the two-point function (0 a 0b) and M. a + M. b = M. c for the three-point 
function (0 a 0{,0*). These are examples of the Bargmann superselection rules |4[] and follow already 
from Galilei invariance 
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]. It follows from the Bargmann superselection rules that no Galilean 
scaling operator can be hermitian unless it is massless. The 'mass' M. therefore plays quite a different 
role in (non-relativistic) Galilean theories as compared to relativistic ones. It no longer measures 
a deviation from criticality, but should rather be considered as the analogue of a conserved charge. 
Finally, the explicit from of the scaling functions in ( p.24 , 2T25|) depends on the way the Galilei 
transformation is realized.^] 

So far, we have always considered both space and time to be infinite in extent. In some appli- 
cations, however, one is interested in situations where the system is 'prepared 'at t — and is then 
allowed to 'evolve' for positive times t > 0. We must then ask which subset of the Schrodinger 
transformations will leave the t — boundary condition invariant as well. Indeed, inspection of the 
generators ( |2.20| ) shows that the line t = is only modified by X_i and that furthermore the subset 



{A 0j i, Y±\j 2 i M } closes. We may therefore impose the covariance conditions (|2.23|) with n = 0, 1 and 



m 



±1/2 only ||53|| . Then the two-point function is 



>dh-r\ )<4(/- 2 .r 2 )) ■ O i2 <\u,..\h [ J 



(x2-xi)/2 



(*i - h)~ Xl exp 



Mi (ri - r 2 ) 

2 ti-t 2 



; t 1 >t 2 >0 
(2.26) 

Compared to eq. fl2.24p , there is no more a constraint on the exponents X\ j2 , because time translation 
invariance was no longer assumed. 

Although Schrodinger invariance of n-point functions was imposed at the beginning of this section 
ad hoc, there exist by now quite a few critical statistical systems with 8 = 2 where the predictions 
( p.24| , pT25| . [2.26D have been reproduced [^, |51f| . Models where some Green's functions coincide with 
the expressions found from Schrodinger invariance include the kinetic Ising model with Glauber dy- 

the symmetric and asymmetric non-exclusion 



namics 



43"! , the symmetric exclusion process |J7], |1 1 



processes |[L01]| , a reaction-diffusion model of a single species with reactions 2 A <-> 20 [f|9| and in 



3 For example, one may modify the generators (2 
lattice constant a 



and (^2 



to emulate the effect of a discrete lattice with 
and conformal 



This works for free fields for both Schrodinger |54[] and conformal [p6| invariance. In the context 
of conformal invariance, the best-known example of the dependence of the correlators on the realization are the 
logarithmic conformal held theories, see B^, |93] for recent reviews. 
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the axial next-nearest neighbour spherical model (ANNNS model | 107|| ) at its Lifshitz point fiO |. In 
section 5, we shall consider in detail the phase ordering kinetics of the 2D and 3D Glauber Ising 
model and several variants of the kinetic spherical model with a non-conserved order parameter, see 
|58|| . Finally, for short-ranged interactions, one may invoke a Ward identity to prove that invariance 
under spatio-temporal translations, Galilei transformations and dilatations with 8 = 2 automatically 



imply invariance under the 'special' Schrodinger transformations [53 



3 Infinitesimal local scale transformations 



We want to construct local space-time transformations which are compatible with the strongly 
anisotropic scaling (O) for a given anisotropy exponent 6^1. The first step in such an under- 



taking must be the construction of the analogues of the projective conformal transformations 
and the Schrodinger transformations ( |2.15| ). This, and the derivation and testing of some simple 
consequences, is the aim of this paper. A brief summary of some aspects of this construction was 
already given in [ER B7|. The question whether these 'projective' transformations can be extended 



towards some larger algebraic structure will be left for future work. 



For simplicity of notation, we shall work in d 
d > 1 will be obvious. 



1 space dimensions throughout. Extensions to 



3.1 Axioms of local scale invariance 

Given the practical success of both conformal (6 = 1) and Schrodinger (9 = 2) invariance, we shall 
try to remain as close as possible to these. Specifically, our attempted construction is based on the 
following requirements. They are the defining axioms of our notion of local scale invariance. 

1. For both conformal and Schrodinger invariance, Mobius transformations play a prominent 
role. We shall thus seek space-time transformations such that the time coordinate undergoes 
a Mobius transformation 

t ^ t , = at + P a5 _p 1 = l (31) 
7t + o 

If we call the infinitesimal generators of these transformations X n , (n = —1,0, 1), we require 
that even after the transformations on the spatial coordinates r are included, the commutation 
relations 

[X n , X m ] = (n - m)X n+m (3.2) 

remain valid. Scaling operators which transform covariantly under ( |3.1|) are called quasipri- 
mary, by analogy with the notion of conformal quasiprimary operators ||. 

2. The generator Xq of scale transformations is 

X = -tdt - - Q rd r - ~ (3.3) 
where x is the scaling dimension of the quasiprimary operator on which Xo is supposed to act. 

3. Spatial translation invariance is required. 

4. When acting on a quasiprimary operator <f), extra terms coming from the scaling dimension of 
(f) must be present in the generators and be compatible with ( |3.3| ). 
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5. By analogy with the 'mass' terms contained in the generators ( |2.20|) for 6 = 2, mass terms 
constructed such as to be compatible with 9 ^ 1,2 should be expected to be present. 

6. We shall test the notion of local scale invariance by calculating two-point functions of quasipri- 
mary operators and comparing them with explicit model results (see section 5). We require that 
the generators when applied to a quasiprimary two-point function will yield a finite number of 
independent conditions. 

The simplest way to satisfy this is the requirement that the generators applied to a two-point 
function provide a realization of a finite-dimensional Lie algebra. However, more general ways 
of finding non-trivial two-point functions are possible. 



3.2 Construction of the infinitesimal generators 

The generators X n which realize ( |3.2| ) will be of the form 

X n = X<P+XM+Xg") (3.4) 



where Xn = —t n+1 dt is the infinitesimal form of fl3.1|) , Xn 1 ^ contains the action on r and the scaling 
dimensions while X n in ' will contain the mass terms. From Xn \ we already have the commutation 
relations (|3.2|) and Xn' Iir> will be constructed such as to keep these intact. 

We now find X n U) . Since for time translations, X_\ = —d t and using (|3.3|), we make the ansatz 



X n = -t n+1 d t - a n (t, r)d r - b n (t, r) (3.5) 

and have the initial conditions 

T X 

a_i = 6_i = , a (t,r) = - , b (t,r) = - (3.6) 

To have consistency with (|3.2j), we set first m = —1, yielding [X n , = (n+ l)X n _i. This gives 
the conditions ^ 

(n + l)a n _i , — ^ = (n + l)6 n _i (3.7) 



dt v ' ' dt 



with the solutions 



ik(r)t n - k 



«n(t,r) = ( H k ^\ J A k 

k=0 ^ ' 

K(t,r) = Y.(lX\)B k {r)r- k (3.8) 

! — n \ / 



where A n (r) and B n (r) are independent of t and Ao(r) = r/6 and Bo(r) = x/9. Next, we set m = 
in (|3.2| ), yielding [X n , X ] = nX n and thus obtain the conditions 

1\ da n 1 da n , db n 1 . , 

n + -e) an = t -dt + 9 r -fr ' n6 " = t ^ + ^ (3 ' 9) 

Inserting ( p.8| ), we easily find 

A fc (r) = A fco r efc+1 , S fc (r) = 5^ (3.10) 
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where A k o, B k0 are constants and Aqq = 1/6, B 00 = x/8. The values of these constants are found 
from the condition [X n , X\] = (n — l)X n+ i. Using the explicit forms a\(t, r) = 28~ 1 tr + A\or e+l and 
bi(t,r) = 2x9~ 1 t + B w r 9 , we obtain 



t I t-^r- + - I r— a n — t r ) + A W U | r- 



dt 9 \ dr 



dr 



+ IK 



, db n 2 ( db n \ db n _ x 

t ( t— + - [r— - xt j + A 10 Ur— - B 10 9Ur a n 



(n - l)a n+ i 
(n - l)b n+1 



(3.11) 



where U = r /t. Using ( |3.9| ), the first of these becomes 
(1 da n ( 1\ 2 TT ( da n 



+ l)fln 



(n — Da. 



n+l 



(3.12) 



Insertion of the explicit form of the a n known from fl3.8| , |3.10| ) leads to (terms with k = 0, 1 cancel) 



E 



k=2 



n + l 



(n + k) 



n + 2 



l, + l J 1 " \ k + l) {n - l) ) Ak » + \ k 
+ [(n - 1) (M n0 A 10 - A n+lj0 )] f/™ +1 = 

which must be valid for all values of U. This leads to the conditions 



til H t l \(k-2)A k _ lfi A w 



(3.13) 



n+l 
k + 



1 Un + k) 



n + 2 
k + 



1 J (n - 1) J A k0 + 
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n + l 
k 



{k - 2)A k _ lfi A 
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(3.14) 



9A nQ A l0 ■ Vn>2 (3.15) 
From ( |3.15| ), we have A n0 = 9 n ~ 2 A 2 oA™q 2 for all n > 2. Inserting this into ( |3.14j ), this is automatically 



satisfied for all k > 2 because of the identity 

(»+*)-(fc+i) (n - 1) + ("* 1 ) (<: - 2) = 



n + l 
jfc + 1 



(3.16) 



In particular, A 20 remains arbitrary. Finally, using the identity 

n + l 



E 

k=2 



k + l 



(l+x) n+1 -l . 1 . 

x k = - (n + l) n(n + l)x 

x 2 



(3.17) 



and using the initial conditions for a n , we obtain the closed form 



I — n \ / 



k=0 

n + l 



t n r + -n(n + l)A w t n ~ v r 



A 



20 



9A 2 



+ 



A 



20 , t n+l r l- 



(9A W 



(3.18 

(1 + 9A 10 r e /t) n+1 - 1 



which depends on the three free parameters A w , A 20 and 9. 
Next, using ( |3.9|) , the second of the relations ( |3.11|) becomes 



+ nb n - ^-t n + A 10 Ur^ - BxoUr-'a^ = (n - l)b n+1 (3.19) 
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which we now analyse. Inserting the form ( |3.8|) for a n and b n , we obtain the condition 

n 

£ 



k=i 



k + i ) (n + k)B k0 + 6 I n + ) ((/, - DD^uAiu - A,,- SM D SI) , 



A . + , j ■ ,-„,, ^ + ((n + l)n - (n- l)(n - 2) - 2) | (3.20) 

+ [9 (nB n0 A 10 - A n0 B 10 ) - (n - l)£„+i, ] C/ n+1 = 

which must be valid for all values of U. Now the term of order 0(U°) vanishes and the term of order 
0(U n+1 ) yields the recurrence 

B n +i,o — 7 (nB n oAio — A n oBio) 

n — 1 

= (nB n0 A 10 - 6 n - 2 A n w 2 A 20 B 10 ) (3.21) 

n — 1 

If we insert this into the above condition for terms of order 0(U k ), k = 1, . . . ,n and use again the 
identity flOU) , we see that all these terms vanish. The final solution for the coefficients B n o is 

B n0 = (n - 1) (9A w ) n - 2 B 20 - {n - 2)9 n ~ 2 A n 10 ~ 3 A 20 B W (3.22) 

and where Bio and -B20 remain free parameters. Using the identities 



k=2 



A/n + lV, . , ((n-2)a;-3)(l + x) n 3 . 1 , . 

E( jfc + l )( fc - 2 )^ = ^ + - + 2(n + l) + -n(n + l)x (3.23) 

k=2 ^ ' 

and the initial conditions for the b n , the closed form for b n (t,r) reads 

Wl , r) = ^l lf + !!(!!±l) f -. r ' fll ,(i-^ 
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+ t „ A m B 20 -2A 20 B la [(n + 1} + (n _ x) (l + MiorV()1 

+ r+V -, 2^ - ^.B,, [x _ (1 + (3 24) 

,n ^20-^10 /-, . /, A d lr\ n 

and depends on the free parameters A w , A 20 , B 10 , B 20 and 9. 

The results obtained so far give the most general form for the generators X n satisfying [X n , X m ] = 
{n—m)X n+m with m = —1, 0, 1 and n G Z and with X given by (|3T3"|). If we were merely interested in 
the subalgebra spanned by {X_i, X ,Xi}, we could simply set A 20 = B 20 = 0, since those parameters 
do not enter anyway in these three generators. 

We now inquire the additional conditions needed for [X n , X m ] — (n — m)X n+m to hold with 
n, m G Z. Given the complexity of the expressions ( p.!8|J3.24D for a n and b n , respectively, it is helpful 
to start with an example. A straightforward calculation shows that 

[X 3 , X 2 ] = X 5 + 9 2 r 59 (9Aj - A 20 ) {A w A 20 rd r + (4A 10 B 20 - 3A 20 B 10 )} (3.25) 

The extra terms on the right must vanish. This leads to the distinction of four cases which are 
collected in the following table 
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^20 


-Bio 


B 2 


1. 


^ 


OAL 

WJ A 10 


^ 


4 


2. 










7^0 


3. 












4. 





7^0 





7^0 



and we now have to see to what extent these necessary conditions are also sufficient. Indeed, straight- 
forward but tedious explicit computation of the commutator [X n ,X m ] shows that it is equal to 
(n — m)X n+rn in all four cases. While for case 1 the expressions for a n and b n are still lengthy, they 
simplify for the three other cases 



_1 (n + l)t n r 
a n (t, r) = { 9-\n + l)t n r + \n(n + l)^ 1 ^ 1 A 10 
9-\n + l)t n r + |(n 3 - n)t n ~ 2 r 29+1 A 20 



case 2 
case 3 
case 4 



b n (t,r) 



e 



n{n+l) , n -l p 1 n A -n + n-2„26 p 
1 r £> 10 H T—t r tf 2 Q 



xt n + 



n(n+l) , n -i h p 
n'—n ^n— 2^.26 Q 



10 



20 



case 2 
case 3 
case 4 



(3.26) 



(3.27) 



We summarize our result as follows. 
Proposition 1: The generators 

X n = -t n+1 d t - a n (t, r)d r - b n (t, r) 



(3.28) 



where 



a n (t,r) 



H+l t n r + ^n(n + l)A 10 f 
A 



n-i r e+i 



1 - 



i20 



9A 2 



20 /" ' 



and 



b n (t,r) = 



n+1 



xt n + 



(^io) 3 

n { n + l ),n-l 



(l + 9A w r e /t) n+1 -l 



(3.29) 



10 



20 



)+«*"^(l + «>W/t)" 



A10-B20 — 2A20-BlO V, M v . / iWl i d/l 0/ + \ n l 

^743^ K n + 1) + (n - 1) (1 + 0A w r /t) \ 

-fl2Aioi?20 — 3A20-B10 



+t" +1 r 



and where one of the following conditions 



[l-(l + 9A w r 6 /t) n ] 



6AL , B 10 ^0 , 5 20 ^0 



(3.30) 



(1.) A 10 ^0 , A 20 -^ 10 

(2.) A 10 = A 20 = , B w ± , £? 20 ^ 

(3.) A w ^ , A 20 = , Bio 7^ , B 20 = 

(4.) A 10 = , A 20 ^ , S 10 = , B 20 + 



(3.31) 



holds, are the most general linear (affine) first-order operators in d t and d r consistent with the axioms 
1 and 2 and which satisfy the commutation relations [X n ,X m ] = (n — m)X n+m for all n,m e Z. // 



12 



only the subalgebra {X±i } is considered, A w and B w remain arbitrary and the generators X± lt0 are 
those of case 3. 

Before we construct the mass terms, we consider space translations, generated by — d r . There 
will be a second set of generators 

Y m = + (3.32) 

where Y^ contains the action on r and Ym 11 ^ contains the mass terms. Given the form Q2.2(J| ) of the 
generators of the Schrodinger algebra, we do not expect any terms proportional to d t to be present 
in the Y m . Indeed, if we tried to include terms of this form, it is easy to see that one were back to 
the case 6 = 1, that is conformal invariance. 

The following notation will be useful. Let 

9 = 2/N (3.33) 
which defines N. We write (up to mass terms to be included later), 

2 (da k (t,r) n db k (t,r)\ 

y - = ^ = -j75rPi)(-5^* + -w J J (3 ' 34) 

where m = — y + k and k is an integer. Here, a n and b n are those of the proposition 1. In particular, 
Y-n/2 = -d r and Y n _ N/2 is obtained from [X n ,Y_ N / 2 ] = \N(n + l)Y_ N/2+n . If A w , A 20 , B w , B 20 
would all vanish, we have indeed [X n , Y m ] = {\Nn — m)Y n+m and we now look for the conditions on 
the parameters which will retain this commutator for all values of n and m. 

For the general situation given by (|3.31| ), direct calculations show that 

[X-x, Y m ) = (-j - rr?j Y m _ x , [X , Y m ] = -mY m (3.35) 

throughout, but the commutator with Xi is more complicated. We shall consider the four cases one 
by one. 

1. For case 1, we consider 

K lk := [X u Y k . N/2 ] -(N-k) Y k+1 ^ N/2 (3.36) 
Since this is still very complex, we expand in A w and find 

K lk = -^^-kB^t^r 26 - 1 - 9 2 kA 10 B 10 t k - 1 r 2e ~ 1 

+ 62{7 ~ 5d) k(k - l)t^V 3e -M 10 5 10 + O (A 2 10 ) (3.37) 

Therefore, if A w , B w , B 20 are all independent, it follows that B w = B 20 = (the other possibility 
A w = B 20 = reduces to a special case of either case 2 or 3 and will be treated below). In this case, 
we expand further and find 

K lk = - e(4H1 6 )( ^ 1) ^m + O (Af ) (3.38) 

Therefore, Aiq = and the case 1 has become trivial, unless 6 = 1. On the other hand, for 6 = 1 
there is a non-trivial solution of the K\ k = 0, namely B 20 = ^A W B 10 . It is now straightforward to 
check that the algebra of the generators X n , Y m indeed closes for all values of n and m. 
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2. For case 2, we have 

K lk = -^^H^-V*- 1 (3.39) 

which implies either B 2 q = for generic 9 or else 9 = 2 and B 2 q arbitrary. The remaining commutators 
[X n , Y m ] are equal to (nN/2 — m)Y n+m for both possibilities. 

3. For case 3, we have 



K lk = --^±^-kA\ t k - x r™d r - e 2 kt k - 1 r 20 - 1 A lo B lo 

which implies A10 = 0. We therefore recover the case 2. 
4. Finally, for the case 4, we have 



(3.40) 



(3.41) 



and therefore for generic 6 1 , we must have v4 20 = B 20 = which is trivial or else we must have 9 = 2. 



In that last case, we consider 



[X 2 , Y k „ N/2 ] = QiV - kj r fe+2 -iv/2 - \k(k - l)A 2 20 t k ~ 2 r 8 d r - h(k - l)A 20 B 20 t k - 2 r 7 (3.42) 

and the extra terms on the right only vanish if A 2Q = 0. This reproduces a special situation of case 
2. 

In conclusion, the unwanted extra terms in [X n , Y m ] are eliminated in three cases, namely 

(i) N generic , B w ^ , B 20 = , A 10 = , A 20 = 

(ii) N = l , Bio 7^ , B 20 ^ , A 10 = , A 20 = (3.43) 
(hi) iV = 2 , S 10 ^ , B 20 = ^A W B 10 , A 10 + , A 20 = ^ 

For the three cases ( |3.43| ) we list the explicit form of the generators X n with n G Z and F m with 
m = k — N/2 and /c e Z in table |I[ In all three cases, the generators depend on two free parameters. 

We still have to consider the commutators [Y m , Ye\. Indeed, in the first case ( p.43|) , the commutator 
between the Y m is non-vanishing 

[Y m , Y t ) = (m-£)(N- 2)^^+Ar-i r 2/iv-2 (3 _ 44) 

Unless N = 2/(2 + n) = 1, |, |, |, . . . or N = 2, that is 9 = 1, 2, 3, 4, . . ., there will be an infinite 

series of further generators. In the second case (|3.43|) , there are three series of new generators Zn \ 
i = 0, 1, 2, see below. Finally, in the third case (|3.43| ), the commutator [Y m , Yg] = A w (m — £)Y m+ g, 
see below. Our results so far can be summarized as follows. 



Proposition 2: The generators X n defined in eq. (3. 2A ) with n G Z and the generators Y m defined 
in eq. ( 3. 3$ with m = —N/2 + k and k G Z and where a n and b n are as in proposition 1 satisfy the 
commutation relations 

[X n , X n >] = (n - n')X n+n , , [X n , Y m ] = (^n^- - rnj Y n+m (3.45) 



14 



Table 1: Generators X n and Y^-nii without mass terms and with n, k G Z according to the conditions 
(ii), (hi) ofeq. (CT). 



(i) 


x n 




= t n+1 d t n fNt n rd r {n+ 2 1)x Nt n n{n + l) B w t n ~ l r 2 l N 




Y k . 


N/2 


= - t kQ r _ ^ kBlot k-l r -l+2/N 


(ii) 


X n 




= t n+1 d t \{n + l)t n rd T \{n + l)xt n n{n + x) B 1Q t n ~ l r 2 {n2 ~ 1)n B 2Q t n ~ 2 r 4 




Y k - 


1/2 


= -t k d r - 2kB 10 t k - 1 r - \k{k - l)B 20 t k ~ 2 r 3 


(iii) 


X n 




= -t n+1 d t - A^[(t + A 10 r) n+1 - t n+1 ]d r - (n + l)xt n A 10 r) n - t n ] 




Y k - 


1 


= -(t + A 10 r) k a r -^B 10 {t + A 10 r) k - 1 



in one of the following three cases: 

(i) B w arbitrary, A w = A 20 = B 20 = and N arbitrary. 

(ii) B\o and B 2 q arbitrary, A±o = A 2 q = and N — 1. In this case, there is a closed Lie algebra 
spanned by the set {X n , Y m , Z n 2 ' , Zm\ Zn^} of generators where nfZ and m G Z + \ and 



Z f := -nt n ~ l r 2 , Z^ : = -2t m - l ' 2 r , := 



-2t r 



(3.46) 



J m+m' 

-n'Z^ , [Y m ,Z^]=-nZ 
' n 



with the following non-vanishing commutators, in addition to 
[Y m , Y m ,\ = (m - m') (±B w Z^ +m , + B 10 Z (0) 

Y 7& 
[Xn, Z^] 



(1) 

n+m J 



+ m ) Z, 



r(l) 

J n+m i 



Y 7 (1) 



-z 



(o) 

ra+m' ' 



V 7 (0) 



-n'Z^ (3.47) 



where n,n' G Z and m,m' G Z + | . TTie Lze algebra structure is determined by the parameter 
B\o/B 20 . 

(iii) A w and B w arbitrary, A 20 = Af , B 2Q = ^A W B W and N = 2. Then for all n,m G Z one has 



\Y n , Y^ 



A w (n - m)Y n 



(3.48) 



[X n , X m ] = (n- m)X n+m , [X n , Y m ] = (n - m)Y n+m , 

The verification of the commutators is straightforward. 

For case (i), if N G N and B w = 0, there is a maximal finite-dimensional subalgebra, namely 
{-X"±i,0j Y_ n / 2 , Y_ n / 2+ x, . . . , Y +N / 2 }. For case (ii), the maximal finite-dimensional subalgebra is span- 
ned by {X±i j0 , Kti/2, 4:B 2Q Z^ + B w Zq}. The Schrodinger algebra eq. (|2.21| ) is recovered for AT = 1, 
-Bio = M./2 and B 20 = 0. The inequivalent realizations of the Schrodinger algebra are classified in 



74| and two distinct realizations were found. The first one of that list j74] is the one discussed here 
and the second realization is excluded by our axiom 1. For N = 2 in case (i), the conformal generators 
will be fully recovered once the mass terms have been included. Finally, case (iii) is isomorphic to 
the conformal algebra fl2.6p through the correspondence X n = £ n + I n , Y m = A 10 £ n . 

We now construct the mass terms contained in xi 111 ^ and Ym H \ For us, a mass term is a 
contribution to the generators which generically is not proportional to a term of either zeroth or first 
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order in d t or d r . The preceeding discussion has shown that the terms merely built from first order 
derivatives dt, d r or without derivatives at all have already been found. The simple example outlined 
in appendix C rather illustrates the need for 'derivatives' 8°: of arbitrary order a. For our limited 
purpose, namely the construction of generators which satisfy eqs. ( |3.45| ), we require the operational 
rules 

d? +b = d?d b , [d?,r] = ad?- 1 (3.49) 

together with the scaling d^f(Xr) = A a <9^ r /(Ar) and that for a = n G N, we recover the usual 
derivative. However, the commutativity of fractional derivatives is not at all trivial and several of the 
existing definitions, such as the Riemann-Liouville or the Griinwald-Letnikov fractional derivatives, 
are not commutative |7^, |9T|, On the other hand, the Gelfand-Shilov [[12], |9T| or Weyl ]79" 



fractional derivatives or a recent definition in the complex plane based on the Fourier transform ||115 



do commute. To make this paper self-contained, we shall present in appendix A a definition which 
gives a precise meaning to the symbol <9" and allows the construction of X n and Y m to proceed. In 
the sequel, the identities (|A5| , |A6| , [A7| , [A9| , |A 1 7| ) will be used frequently. 

For generic N, setting x = and B\q = for the moment, we make the ansatz 

X n = -t n+ % - j (n + l)t n rd r - A n {t,r)8 a t (n) - B n (t,r)d b r ^ - d c T ^C n (t,r) (3.50) 

where the functions A n ,B n ,C n and the constants a{n) , b(n) , c(n) have to be determined. From the 
condition [X n , Xq] = nX n , we find the equations 



td t + yrSr ) A n (t. r) - r/(//).4„(/. /') 



N \ N 
td t + -j rd r) B n(t,r) - -b(n)B n (t,r) 



N 



N 



td t + Y rdr J C '™( t ' r ) ~ -^ c i. n ) c nit,r 

with the solutions, where u = r 2 ' 7V t _1 

A n (t,r)=t n+a ^A n (u) , B n (t,r)=t n+Nb ^ 2 B n (u) , 



nA n (t,r) 
nB n (t,r) 
nC n (t,r) 



C n {t,r) 



t 



n+Nc 



(n)/2 C 



u 



(3.51) 



(3.52) 



f n-l+a(n) ( ( n + a ( n )) ^ 
j.ri-l+Nb(n) /2 
£n~l+Nc{n)/2 



N 

n + —b(n) ) B n 



n + —c(n) ) C r , 



i and find 






r 2/N \ 
f An) 


(n + l)t n ~ 


-l+o(n-l) a , 




(n + l)t n ~ 


-l+JV6(n-l)/2*3 

*^n— 1 


T c ») - 


(n + l)t n ~ 


-l+ATc(n-l)/2/7 

Wi— 1 



(3.53) 



where the prime denotes the derivative with respect to u. Since this must be valid for all values of t 
and r (or t and u), we find 



a(n) = a(n 



b(n) = b(n - 1) = ^ , c(n) = c(n - 1) = ^ 



(3.54) 
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where a, b, c are n-independent constants, and 

(n + a) An(u) — uAn(u) = (n + l)A n -i(u) 
{n + b)B n {u)-uB n '(u) = (n + l)B n _i(u) 

(n + c)C n (u) -uC n '{u) = (n+l)C n _i(u) (3.55) 
In addition, we have the initial conditions 

A-i(u) = Ao(u) = B- X (u) = Bq(u) = C_i(w) = C (u) = (3.56) 
The solution of this is, e.g. for A n , 

xw = E«4 HI K +a ( 3 - 57 ) 
fc=i ^ ' 

where the are free parameters. Similar expressions hold for B n and C n , where a is replaced by b 
and c, respectively, and free parameters fa, Ik are introduced. 

In the sequel, we shall concentrate on quasiprimary operators which are assumed to transform 
covariantly under the action of X±i only. We repeat the explicit expression for the generator X\ of 
'special' transformations, in the simplest case, 

Xx = -t 2 d t - Ntrd r - ar 2 ^ N d a t - fa 2 ^ N d 2b / N - ld ^/N r 2(i+c)/N {3 58) 

and where a, fa 7 are free parameters. 

For the physical applications, it is now important to check the consistency with the invariance 
under spatial translations, generated by Y_ N / 2 = —d r . In particular, from eq. ( |3.45| ), we should have 
[Xi, Y_jv/2] = NY^ N /2+i- From this, we easily find 



Y_ N/2+l = -td r - ^(1 + ay^^d? - + h)r ^)/N-l d 2 b /N _ |7 (1 + c)d 2c/N r 2(l +C ),N-l 

i\ z N z N z 

(3.59) 

Acting again on this with Y_^/2, we have the commutator 



[Y-n/2 + i, Y„ N/2 ] = -^(1 + a)(l + a - N/2)r 2 ^' N - 2 dt 

_i^(l + 6 )(1 + b _ N/2 y(l + b)/N-2 d 2 b /N _ 47 (1 + c)(1 + c _ N/2)g ? c /N r 2(l +C )/N-2 

and a sequence of further generators may be constructed through the repeated action of Y_ N / 2 . The 
number of these generators will be finite only if the conditions 

^(1 + a) = h e N , ^(1 + b) = k 2 e N , ^(1 + c)=k 3 eN (3.60) 

are satisfied. That means that the realizations under construction will be characterized by the value 
of N and the three positive integers k\, k 2 , k%. We shall call the ki the degrees of the realization. 

A further consistency check, for N integer, is provided by the condition [Xi, Y N / 2 ] = or equiva- 
lently, for the N + 1-th iterated commutator 

[X 1 [---[X 1 ,Y_ N/2 ]---]] =0 (3.61) 
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Direct, but tedious calculations show that this is satisfied if either (i) a ^ and (3 = 7 = or 
alternatively (ii) a = and 7 7^ 0. We call the first case Typ / and the second case Typ II. 

If we consider the generators of Typ I, we see that for N — 1 and k% = 2, we recover the generators 
( p.20| ) of the Schrodinger algebra, with a = Ai/2. This explains the origin of the name 'mass term' 
for the contributions to X n , Y m parametrized by a, (3, 7. Furthermore, for N = 2 and k\ = 2, let 

z = t + \far , z = t — yfax (3.62) 

with a = — 1/c 2 where c is the 'speed of light' (or 'speed of sound'). Therefore X n = £ n + £ n and 
Y n = i(£ n — £ n ) where the conformal generators £ n ,£ n are given in (|2.8|). Usually one sets c = 1 and 
the presence of a dimensionful constant is then no longer visible. The Schrodinger algebra generators 
( p.20| ) are also recovered for Typ II with /c 2 = k 3 = 2 and (3 + 7 = M.J2. These special cases already 
suggest that the choice ki = 2 may be particularly relevant for physical applications. We find a third 
example with a finite-dimensional closed Lie algebra in the presence of mass terms for Typ II with 
N = 2 and degree k^ = k 3 = 2. Then the commutators read 

[X n , X m ] = (n- m)X n+m , [X n , Y m ) = (n - m)Y n+m , [Y n , Y m ] = (p + j)(n- m)Y n+m (3.63) 

and we recover case (iii) of the proposition 2, after identifying A 10 = (3 + 7 and B w = 27. 
From now on and for the rest of this paper, we shall always take B w = 0. 



3.3 Dynamical symmetry 



For k\ = &2 — k s = 2, our realization acts as a dynamical symmetry on certain linear (integro-) diffe- 
rential equations with constant coefficients as we now show. In d spatial dimensions, generalizing the 
above constructions along the same lines as in appendix B, we consider the generalized Schrodinger 
operator 

y 



s 



-ad? +[~—\ d r d r 



(3.64) 



and the generators X_i = —d t , Y 



(0 

N/2 



-d ri together with the Typ I generators with k\ 



Y 



X 

(i) 

■N/2+1 



o N n N 

-ta t r ■ o r X 

2 2 

-t 2 d t - Ntr ■ d r - Nxt - ar 2 d t 
N 



2aN-l 



(3.65) 



with i = 1, . . . , d and we have 



fax- 



's, Y (l) 



N/2 



S, Y (t) 



N/2+1 



(3.66) 



which shows that 5 is a Casimir operator of the 'Galilei'-type sub-algebra generated from 



{X-i, Y-h/2+ 1} as gi ven i n ( P-65|) . Furthermore, 



[s,x 



-NS , [S, Xi] = -2NtS + aN 2 [x - - + 



d N — 1 



N 



(3.67) 
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In addition, since [X%, X_jv/ 2+ J = (N — k)Y^, 2+k+1 , it follows immediately that [S,Y_ 
all k 7^ N, since from the Jacobi identities 



'(0 



(0 l _ n f or 

AT/2+feJ — u 1UI 



(AT — A;) -1 



> 1 -N/2+k 



Y 



N/2+k' 



[S,X 1 







(3.68) 



on the solutions of <S"0 = and [<S,Xi]-?/> = by induction over k. Additional generators created 
from the commutators [Y m , Y m /] are treated similarly. Therefore, we have shown: 



Proposition 3: The realization ( P-6d{) of Typ I generated from {X-x, Xi, Y^, 2 }, i = 1, . . . , d sends 
any solution ip(t, r) with scaling dimension 



d N-l 



of the differential equation 



X= 2 



N 



(3.69) 



N 



-ad? +[- ) d r d r )V(t,r) = 



(3.70) 



into another solution of the same equation. If we construct a free-field theory such that (|3.70|) is 
the equation of motion, then x as given in ( |3.69j ) is the scaling dimension of that free field ip. That 
theory is non-local when A" is not a positive integer. 

Similarly, for Typ II with fc 2 = ks = 2 and d — 1 for simplicity (we shall refer to this case in the 
sequel as Typ Ila), we consider 

S = -0J + 7)a + ifif (3.71) 



and the generators of (|3.65| ) are replaced by (k 2 = k 3 = 2) 



1 x 
X = -td t - -rd r - - 

X 1 = -t 2 d t - -trd r -±t-{P + 7 )r 2 9 r 2 ~ e - 7 2(2 - - 7 (2 - 6){l - 6)d; B (3.72) 

a o 



3I-6 



y_ w/2+1 = -td r -(l3 + ^Ord 2 - -^{2-6)8, 

Again, ( |3.66| ) holds so that S is a Casimir operator of the 'Galilei' sub-algebra generated from 
{X-x, Y_ N /2, ^-jv/2+i}- In addition, instead of (|3.67| ) we have 



[S, X ] = -S , [S, X,} = -2tS + £±I f 2x - + 1 - (2 - 0)) (3.73) 

Therefore, we have the following dynamical symmetry. 

Proposition 4: The realization ( 3. 7£\ ) of Typ Ila generated from {X-i,Xi, Y_^/ 2 } sends any solution 

6-1 



ij}(t,r) with scaling dimension 
of the differential equation 



x 



+ 



7 



(3 + 1 



(2 



S^(t, r )=[-(j3 + 7)$ + -d° r ) r) = 



(3.74) 



(3.75) 
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into another solution of the same equation. This means that the ratio (3/^ is a universal number and 
will be independent of the irrelevant details (in the renormalization group sense) of a given model. 
As before, x is the scaling dimension of the free field whose equation of motion is given by ( |3.75| ). 

While these dynamical symmetries were found for k\ = k 2 — ^3 — 2, there is one more possibility 
for Typ II with k 2 = k 3 = 3 (to be called Typ lib in the sequel). Consider 

S = -3(3 - 0) 7 $ + j 2 d e r (3.76) 
and the generators now read (A; 2 = /c 3 = 3) 

1 x 
X Q = -td t - -rd r - - 

X x = -t 2 d t -ltrd r -^t-(3r 3 d 3 r - e -79 r 3 -V (3.77) 


If we take /3 + 7 = 0, we recover indeed eqs. (|3.66|) so that S is again Casimir operator of the 'Galilei' 
sub-algebra and 

[S, X ) = -S , [S, X 1 ] = -2tS + ^^7 (x - (3.78) 
and we have the following statement. 



Proposition 5: The realization ( 3. 77\ ) of Typ lib with k 2 = k% = 3 and (3 = —7 sends any solution 



ip(t,r) with scaling dimension x = 1/2 of the differential equation 

Sip{t, r) = (-3(3 - 0) 7 $ + 1^ r) = (3.79) 

to another solution of the same equation. It follows that for Typ II there are two distinct ways of 
realizing a dynamical symmetry, if 7^ 2. 

All possibilities to obtain linear wave equations with constant coefficients from Casimir operators 
of the above simple form are now exhausted. We illustrate this for Typ I. A convenient generalized 
Schrodinger operator S should satisfy [S, X-i] = [S, Y^n/ 2 ] = and this implies S = S(d t , d r ). Using 
the identity (|A18 ) and writing S = S(u,v), one has for k\ = k 



& Y -m + ^] = ~Q^d r - - ^ ( e r (3.80) 

1=1 v 7 

which for k > 2 contains several terms with powers of r which cannot be made to disappear and 
therefore [S, Y-n/ 2+ i] = is impossible. Similarly, for Typ II only for the case k 2 = k% = 3 a 
compensation of one more term is feasible. For the other cases, Casimir operators can of course be 
constructed in a straightforward way, but we shall not perform this here. 

For N — 1 (or = 2), Typ I (h = 2), Typ Ha (k 2 = k 3 = 2) and Typ lib (k 2 = k 3 = 3) 
coincide and we recover as expected the known dynamical symmetry of the free Schrodinger equation 



84] , |50| , g. For N = 2 and k\ = 2, Typ I gives the Klein-Gordon equation with its dynamical 
conformal symmetry. Finally, for iV = 2, k 2 = k% = 2 and f3 + 7 = A 10, Typ Ha is identical to the 
generators (|3.28| , |3.34] ) with the identification B w = 27. 
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From the different forms of the wave equations ( |3.70| , |3~7r5|) we see that Typ I and Typ II describe 
physically distinct systems. The propagator resulting from Typ I, which in energy-momentum space 
is of the form (u N + p 2 ) -1 , is typical for equilibrium systems with so strongly anisotropic inter- 
actions, that the quadratic term ~ uo 2 which is usally present is cancelled and the next-to-leading 
term becomes important. This occurs in fact at the Lifshitz point of spin systems with competing 
interactions, as we shall see in a later section. On the other hand, the propagator from Typ II, of 
the form (u + p d )~ l in energy-momentum space is reminiscent of a Langevin equation describing the 
time evolution of a non-equilibrium system (furthermore, if we had tried to describe such a real time 
evolution in terms of the propagators found for Typ I, we would have encountered immediate prob- 
lems with causality.) Indeed, we shall see that aspects of aging phenomena in simple ferromagnets 
can be understood this way. 



3.4 The two-point function 

The main objective of this paper is the calculation of two-point functions 

G = G(t u t 2 - r l5 r 2 ) = (Mh, ri)Mh, r 2 )> (3.81) 

of scaling operators <pi from its covariance properties under local scale transformations. We shall 
assume that spatio-temporal translation invariance holds and therefore 

G = G(t,r) , t = t 1 -t 2 , r = n-r 2 (3.82) 

Since for scaling operators invariant under spatial rotations, the two points can always be brought 
to lie on a given line, the case d = 1 is enough to find the functional form of the scaling function 
present in G. To do so, we have to express the action of the generators Xo,i and Y m on G. Each 
scaling operator 0, is characterized by either the pair (x^ctj) for Typ I or the triplet (xj,/?j,7j) for 
Typ II.Q By definition, two-point functions formed from quasiprimary scaling operators satisfy the 
covariance conditions 

X n G = Y m G = (3.83) 

Since all generators can be obtained from commutators of the three generators X±i, K_jv/2> explicit 
consideration of a subset is sufficient. For the three cases considered above, the results are as follows. 

1. For Typ I, the single condition 

a 2 = {-l)- N oc x (3.84) 

is sufficient to guarantee that (|3.83|) is satisfied, together with the covariance under all commuta- 
tors which can be constructed from the X n , Y m . We merely need to satisfy explicitly the following 
conditions 

X G(t, r) = f-td t - y rd r - Nx^j G(t, r) = 

X!G(t, r) = (-t 2 d t - Ntrd r - 2N Xl t - ai r 2 d? _1 ) G(t, r) 

+2t 2 X G(t, r) + Nr 2 Y_ N/2+1 G(t, r) = (3.85) 



Y_ N/2+1 G(t,r) = f-tdr-^rd^G^r) 



4 T he indices oti,f3i,ji refer here to the two scaling operators and have nothing to do with the indices used in eq. 

(EH- 
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where 2x = X\ + x 2 . This makes it clear that time and space translation invariance are implemented. 
If we multiply the first of eqs. ( |3.85| ) by —t and add it to the second one and then multiply the third 
of eqs. (|3.85|) by —2/Nr and also add it, the condition XiG(t,r) = simplifies to 



-Nt (xx -x 2 )G(t,r) = 

which implies the constraint 

Xi = x 2 

The two remaining eqs. ( |3.85| ) may be solved by the ansatz 

G(t,r)=5 XUX2 r-^Q(tr- 2 / N ) 
which leads to an equation for the scaling function f2(w), where v = tr~ 2 l N 

(ati0f -1 - v 2 d v - Nx x ) Q(v) = 
and the boundary conditions (see section 2) 

fi(0) = tt , n(v) ~ n^y-^ 1 for v -> oo 



(3.86) 
(3.87) 

(3.88) 
(3.89) 
(3.90) 



where f2 0iO o are constants. Eqs. ( p.88p.89| , |3.90|) together with the constraints (|3.84| , |3.87|) determine 
the two-point function and its scaling function Q(v) and constitute the main result of this section 
for the realizations of Typ I. 

2. Similarly, for Typ Ha (k 2 = k% = 2), the conditions 



(3 2 = (-1)-^% , 72 = (-1) 



-l- 



7i 



(3.91) 



are enough to guarantee that the quasiprimarity conditions ( |3.83| ) hold and 



X G(t,r) 
X x G{t,r) 



-td t - -rd r - — 



G(t,r) = 

{-t 2 d t - ~trd r - - (Pi + n)r 2 d 2 - e - 2(2 - O^rd 1 -^ G(t, r) 



+2t 2 X G(t,r) + -r 2 Y_ N/2+1 G(t,r) = 



(3.92) 



Y- N/2+1 = (-tdr-OiPx + ^rd 2 - 6 -29(2-9) ll d 1 r - e )G(t,r) = 
In the same way as before, the condition X\G(t, r) = can be simplified into 



-t (xi - x 2 ) G(t,r) = 

u 



(3.93) 



which implies again the constraint ( [3.87 ) . The two remaining eqs. (|3.92|) can be solved by the ansatz 

G(t,r) = 5 XuX2 r 2 ^%(rr 1 / e ) (3.94) 
and lead to the following equation for the scaling function $(w), where u = rt" 1 ^ 9 

(d u + 0(/3i + ^udl" + 29(2 - 9)^01-°) = (3.95) 
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with the boundary conditions 

$(0) = $ , ~ $oowr 2:Cl for u — ► oo (3.96) 

where <3>o, $00 are constants. Eqs. ( |3.94| , |3.95| , |3.96p together with the constraints ( |3.91|J3.87| ) determine 
the two-point function and its scaling function and constitute the main result of this section 
for the realizations of Typ Ha. 

3. For Typ lib (k 2 = k 3 = 3), there are two possibilities. First, we might take 72 = (— l)~ 1+e 7i, 
which would be the same as in eq. Q3.91]) . It turns out, however, that this leads to the same equation 
for the scaling function as for Typ Ha, upon identification of parameters. We therefore examine the 
second possibility 

72 = (-ir 2+ Si (3.97) 
Then, in contrast to the previous cases, we need the additional generator 

M:= Ye t r - JV / 2 + 1 ' Y ~ N ^ = + l)rd ' r ° e ^%~ 9 (3 ' 98) 

We let Pi t 2 = —71,2 and find 
X G(t,r) = (-t8 t -±rd r -j?)G(t,r)=0 

: ' :y '-"' jidfr e r 3 - 2(3 - 9){2 - 9){l - 9) lx d;° ) G(t,r 



( 2 2ri 

X 1 G(t,r) = l-t 2 d t --trd r --^t + ll r 3 d; 

+ \2t 2 X + ^r 2 Y„ N/2+1 - 3r 2 M^) G{t,r) = 
F_iv/2+i = (-tdr - 39(3 - 6) lir d 2 r - e ) G{t, r) + 36r 2 MG{t, r) = 



(3.99) 



In addition to the usual conditions ( ft.83|) for quasiprimarity, we also need explicitly that MG(t, r) = 
0. As we did before, the condition XxG(t,r) = can be simplified and we find 

td t + ±rd r + ^±^jG(t,r) = 

(7i0(2 - 0)(3 - 0) (3rd r + 2(1 - 9)) d; e + fa - x 2 )t) G{t,r) = (3.100) 

(td r + 7i3^(3 - 9)rd 2 r e ) G(t, r) = 

To analyse these further, let H(t,r) := d~ e G(t,r). Then, the last two of the above equations take 
the form, after having also acted with d r on the second equation ( |3.100| ) 



(^39(3 - 9)(2 - 9) (rd 2 r + ~(5 - 29)dA + fa - x 2 )td 1 r +e ^j H(t, r) = 

(tdl +e + -f 1 39(3-9)rd 2 r )H(t,r) = (3.101) 

The only apparent way to make these two equations compabtible is to make one of them trivial or 
else to make them coincide. The first one is possible if 9 = 2 and x± = x 2 and the second possibility 
occurs if 

9 = ^ , x 2 = Xl + l/2 (3.102) 

and these conditions make Typ lib a very restricted one. In order to have X\ 7^ x 2 in either 
Schrodinger or conformal invariance, time translation invariance must be broken |xj but remarkably 
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Table 2: Some properties of three basic types of local scale invariance. The upper part collects 
operator properties and the lower part specific properties of the two-point function. Here S is the 
generalized Schrodinger operator and x is the scaling dimension of the field ip wich solves Sip = 0, 
in d = 1 spatial dimensions. The scaling function for the two-point function satisfies the differential 
equation T>Q(v) = or T>$(m) = 0, respectively. In the boundary conditions, the first line refers to 
u — or v — 0, and the second line to u — > oo or v — > oo, respectively. 





Typ I 


Typ Ha Typ lib 


degree 
masses 

c 

o 

X 


K\ — 1 

/3 = 7 = 

n xN , 1 Aj2f)2 
—C(O t + jl\ o r 

\-{N-l)/N 


k 2 = k 3 = 2 
a = 

i(0-l) + (2-0)/(l + /ty 7 ) 


&2 = &3 = 3 

a = /5 + 7 = 

-3(3 - eyidt + r 2 <9* 
i 

2 


constraints 


X2 = X\ 

a 2 = (—l)~ N ai 


X 2 = Xi 

fa = 72 = -(-l) 9 7i 


X 2 = Xi + |, = | 

72 = ("l) e 7i 


scaling 


G = r~ 2x ^(v) 
v = tr- 2 ' N 


G = t-^ +x ^l°^(u) 
u = r r 1/e 


V 


aid^' 1 — v 2 d v — Nxiv 


d u + 0(A + 7i«" y + 26(2 - 0) 7 i^- y 


a u + 30(3-0)7 lM a 2 - y 


auxiliary 
condition 






9 2 - e $( M ) = 


boundary 
conditions 


n(o) = n 


$(0) = $ 



enough for Typ lib, in spite of the presence of time translation invariance, the scaling dimensions 



of the two quasiprimary operators are not the same. If the above conditions (|3.102|) hold, one can 
make the ansatz 

G{t, r) = 6 Xl+1/2>X2 r^+**V e <S> (rr l / e ) (3.103) 
where satisfies the equation 

(8 U + 7i3fl(3 - d)ud 2 ~ e ) = (3.104) 

with the usual boundary conditions. In addition, we still have to take into account the condition 
MG(t,r) = -(3 - e)2-f l t- ( - xl+X2+2 ~ 9 ^ e d 2 u - 6 ^(u) = 0. Therefore 9 2 " $(m) = and from eq. (EUOll) 



it follows d u <&(u) = 0. The fact that $(w) = const., together with the restrictiveness of this case 
suggests that Typ I and Typ Ha, where fcj = 2, might be more relevant for physical applications. 

In summary, starting from certain axioms of local scale invariance, we have constructed infinitesi- 
mal local scaling transformations. We have shown that these act as dynamical symmetries of certain 
linear equations of motion of free fields. We have also derived the equations which the scaling func- 
tions of quasiprimary two-point functions must satisfy. For easy reference, the main results are 
collected in table |2|. 

Generalizing from Schrodinger invariance, it may be useful to introduce a conjugate scaling op- 
erator <p* characterized by 

(x,(-l) N a) , for Typ I . , 

(*, -(-!)-*/?, -(-1) A) , for Typ Ha ^ Wb) 
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The two-point functions then become 



f S S t~ 2xi Q (tr~ 2/N ) for Tvn I 

at, r) = (m*u r, m „n)) = { f;;;±2lJX^ (ii-v.) ; Z III L < 3106 > 

Only for Typ I with N even, the distinction between (ft and (ft* is unneccessary. For Typ II, we shall 
in section 5 identify (ft* with the response operator (ft in the context of dynamical scaling. 

In general, the infinitesimal generators X n , Y m do not close into a finite-dimensional Lie algebra. 
However, if we apply them to certain states (which we characterized for two-body operators by the 
conditions 03.84| ) and ( |3.91| ) for Typ I and Ha, respectively) all generators built from [Y m , Y m /] vanish 



on these states. One might call such a structure a weak Lie algebra. For the cases considered here, 
it is generated from the minimal set {X_i, X\, Y_ N / 2 }. 

Technically, we have been led to consider derivatives d a of arbitrary real order a. We emphasize 
that all results in this section only use the abstract properties of commutativity and scaling of these 
fractional derivatives, as stated in ( 3.49| ) and formulated precisely in appendix A, see eqs. (|A3,A^, 



They are therefore independent of the precise form ( |A2| ) which is used in appendix A to prove them 
and any other fractional derivative which satisfies these three properties could have been used instead. 
On the other hand, the differential equations ( |3.89| , |3l35D for the two-point functions can only be solved 
explicitly if a particular choice for the action of d a on functions is made. 



4 Determination of the scaling functions 



We now derive the solutions of the differential equations ( |3.89| , |3.95| ), together with the boundary 
conditions (|3.90|J3.9^ ), see also table [| From now on, we make use of the specific definition eq. (|A2|) 
for the fractional derivatives. The comparison with concrete model results will be presented in the 
next section. 



4.1 Scaling function for Typ II 

Having seen that for the realizations of Typ II considered in section 3 only the one of degree k 2 = 
k% = 2 will lead to a non-trivial scaling function, we shall concentrate on this case, called Typ Ha in 
table |2|. We now discuss the solutions of eq. ( |3.95| ), which we repeat here 



(d u + B{Px + li)udl~ e + 29{2 - d)^- 6 ) = (4.1) 

Before we discuss the general solution of eq. ( [4.1[ ), it may be useful to consider the special case 
9 = 1 first. Then the scaling function satisfies 

{d u + (Pi + li)ud u + 2 7l ) $(«) = (4.2) 

The solution is promptly found 

- / $ o (i + (A + 7 i)«r 2/(1+/3l/7l) ; 7i ^ -A iA 3] 

9 W-\ <Mx P (-2 7l u) ; 7i = -A 1 j 

where $o is a constant. Although 9 = 1, this situation does not correspond to the case of conformal 
invariance, as a comparison with the conformal two-point function $ C onf(^) ~ (1 + u 2 )~ x from 
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eq. ( |2.12| ) shows. Furthermore, although eq. ( |3.95| ) does not contain the scaling dimensions explicitly, 
the ratio fti/ji 7^ —1 is indeed a universal number and related to the scaling dimension x\ via 



l + A/71 

as the consideration of the boundary condition eq. fl2.2|) shows. Then we finally have 



$ (1 + 1 U) 

$0 ex P (— 27im) 



7i 7^ -A 
7i = -Pi 



(4.4) 



(4.5) 



We remark that the case /3\ + 71 = is one of the rare instances where Cardy's [18[ prediction of a 
simple exponential scaling function is indeed satisfied. The special case 6 = 1 may be of relevance in 
spin systems with their own dynamics (which is not generated from a heat bath through the master 
equation). This situation occurs for example in the easy-plane Heisenberg ferromagnet, where the 
dynamical exponent z = 1.00(4) in 2D PSfl . The dynamical exponent z = 1 also occurs in a recently 
introduced model of a fluctuating interface |Q and related to conformal invariance. Another example 
with a dynamical exponent z ~ 1 is provided by the phase-ordering kinetics of binary alloys in a 
gravitational field [|9[ . 

We now turn towards the general case. In order to find the solution of eq. (|3.95|) by standard 



series expansion methods, see |79|. |9"]J, we set 



P 
Q 



where p, q are coprime integers and make the ansatz (with cq 7^ 0) 



£ 

n=0 



n/q+s 



(4.6) 



(4.7) 



where the c n and s are to be determined. Although this procedure certainly will only give a solution 
for rational values of 6, we can via analytic continuation formulate an educated guess for the scaling 
function for arbitrary values of 9. Furthermore, since we are only interested in the situations where 
the scaling variable u is positive, we can discard the singular terms which might be generated by 
applying the definition ( |A2|) of the fractional derivative d^. At the end, we must consider the right- 
hand limit u —>■ 0+ and check that indeed lim u ^o+ *&(u) = $0 exists. f] 

Insertion of the above ansatz into the differential equation leads to 



u 



1 £ (h 



v n=0 
00 



s c n u 



n/g 



n=p 



9T [ ^ + s + 1 



n 



(A +7i) ( - + 3-11+2(2-0)7! 



u 



n/q 



and must be valid for all positive values of u. This leads to the conditions s = 0, c n 
n = 1, 2, ... ,p — 1 and, if we set $^ := c p e, to the recurrence 



7l r(g(f-l) + l) 



[M + B] $ 



e-i 



(4.8) 
for 

(4.9) 



5 This formal calculation works here in the same way as for ordinary derivatives since the highest derivative in ( |4.l| ) 
is indeed of integer order. 



26 



where we have set 

A := 0(1+ ft/71) , B:= 3 -29 -ft/71 (4.10) 

To solve this, we set 

= P(£)rV, ; t > 1 (4.11) 
where the function P(ri) is defined by P(l) ■— 1 and for n > 2 by 

n-l 



n r(w ° r(a) (412) 



We first discuss the case ft + 71 7^ 0. Then A ^ and we find the following recurrence for the ai 

a l+l = - lx QA (£ + 1 + B/A) a t ; £>1 (4.13) 
which is solved straightforwardly and leads to 

$t = ( _ 7lM/ r« + i + f>M) 

v ' ; r(^ + i)r(i + 5/A) v ; 

Similarly, for the other case ft + 71 we have A = and we find in an analogous way 

= (- 7l g£)* 1 $ (4.15) 
Therefore the sought-after series solution finally is 

$o^,a(-^ 2 (/9i+7iK) ; 7i 7^ -A 



where 



and the function (E a ,6(z) and the Mittag-Leffier function E a ^(z) are defined as 

hi — k — 



From this is it obvious that the scaling function $(w) as given in ( 4.16 ) does satisfy the boundary 
condition $(0) = $0 — const and has an infinite radius of convergence for 9 > 1 if ft + 71 7^ and 
for 9 > if ft + 7x = 0. Using the identities 



<B hb (z) = (1 - z)- b , £ 2jb ( z ) = 1 F 1 (^b;~-,^j , (£ 2il/2 (^) = e^ 4 



(4.19) 



E 1A (z) = e z , E 1/2;1 (-z) = e z (1 - erf (z)) (4.20) 

it is easily checked that the known solutions ( 4~5l) and (|2.24| ) are recovered, for 9=1 and 9 = 2, 
respectively. 

In figures [I] and ^ the behaviour of the scaling function is illustrated for several values of 9 
and of A. The case ft + 71 = 1 is governed by the properties of the function <£ ajb (z). In figure [l] for 
two fixed values of 9 the effect of varying A is displayed. Since the universal ratio ft/71 is arbitrary, 
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the parameter A as defined in (|4.17|) can take any positive value. Only for 9 = 2 we have A = 1/2 



fixed. The dependence of the scaling function for (3\ + 71 7^ on 9 is shown in figure Qa, where 
we fix A = 1/2. Finally, the scaling function found in the peculiar case (3\ + 71 = is illustrated in 
figure 0b. This case is governed by the Mittag-Leffier function Egi(z) whose properties are reviewed 



in some detail in 91 



In the following cases, the asymptotic behaviour of the scaling function is algebraic, as follows 
from the identities, see |]113|] 



E a , b (-x) ~ — - — -x 1 + 0(x- 2 ) , 0<a<2 
1 (0 — a) 

Cajbi-x) * , rn 1 . x- 1 + 7/ ~ b l ^ b + O (x- 2 , x- 1 -") , Ka<3 (4.21) 
ol (1 — a) 1 (1 — ab) 

In these cases, the physically required boundary condition is satisfied and one may deduce the relation 
between A and the scaling dimension Xi, thereby generalizing (|4.4|) . 



We also see from figure || that on a qualitative level, the cases /?i + 71 7^ and Pi + ji — are 
broadly similar. For sufficiently small values of 9, the scaling function descreases monotonically 
towards zero when u increases. With increasing 9, the scaling function decays faster for large u and 
at a certain value of 9 (at 9 — 2 or 9 = 1, respectively), the asymptotic amplitude vanishes and the 
decay becomes exponential. If we now increase 9 slightly, the scaling function $(w) starts to oscillate. 
Oscillations also arise if for 9 fixed the parameter A is made sufficiently large, as can be seen from 
figure [1]. 



4.2 Scaling function for Typ I 

We write eq. ( |3.89| ) in the form 

(aidf- 1 - v 2 d v - v() Q(v) = 



2xi 

C = N Xl = -f 



(4.22) 
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Figure 2: Scaling functions $(w) for Typ II as given by eq. ( |4.16|) . (a) Several values of 9 as indicated, 
A = 1/2 and (5\ + 71 = 1. (b) Several values of 9 and 71 = —f3\ = 1. 



It is useful to begin with the special case when N is an integer ||55|| . Then the anisotropy exponent 9 = 
2/N = 2, 1, |, |, |, ■ ■ ■ and one merely has a finite number of generators Y m , m = —N/2, . . . , N/2. 
For N = 1 and JV = 2 we recover the scaling functions found from Schrodinger and conformal 
invariance (see section 2) and now concentrate on the new situations N > 3. For N = 4, some explicit 
solutions for a few integer values of ( are given in table [| Given the boundary condition f2(0) = 1, 
these still depend on two free parameters /?, If (3' 7^ 0, these solutions diverge exponentially fast 
as v — » 00 but if we take f3' = 0, we find Q(v) ~ in agreement with the required boundary 
condition, see table |[ 

Using these examples as a guide, we now study the more general case with integer iV and ( 
arbitrary The general solution of eq. ( 3.89j ) for integer iV > 2 is readily found 



n(«) = jr b p v^ p ; T v = 2-^tv-i ( C -p-, 1; 1 + ^, 1 + • • • , ) (4.23) 



p=0 



N ' ' iV' N N ' N N ~ 2 ai 



where 2-PV-1 is a generalized hypergeometric function and the 6 P are free parameters. To be physically 
acceptable, the boundary condition ( |3.90|) must be satisfied. The leading asymptotic behaviour of 



the J-'p for v —>■ 00 can be found from the general theorems of Wright ||113|| (see [[RJ for a brief 
summary) and the asymptotics of Q(v) is given by 



Q ( v ) ~ \ TT^ /_ ,m//v ex P ~ m<n-k v 




N/{N-2) 



which grows exponentially as v — > 00 if N > 2. Clearly, this leading term must vanish, which imposes 
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Table 3: Some solutions Q(v) of eq. ([4.22 ) for N = 4 with O(0) = 1 and their leading asymptotics 
for P' — as v — > oo. Here l u , K u are modified Bessel functions, L u is a modified Struve function, 
erf is the error function (see [Q) and (5,(5' are constants. The abbreviation y := t> 2 /(2 v /a7) is used 
throughout. 



c 



asymptotics 
(for (5' = 0) 



T(3/4) 2 



Vy^U (D+/M1/4 (i)k 1/4 (|) 



/V^r(J) (|) 1/4 [L_ a/4 (y) - (y)] 



+r(I) (I) 1/4 [L-^Cv)-! 



1/4 



e" y + P (4f ) V4 [ e " erf (Vv) - ie "" erf (i VI/) - 2 sinh y] 



l + V^r(|) (1) 3/4 [L 1/4 (y) - I-1/4 
3/4 [l_ 1/4 (y) - V 



+/M/4 (!) 



+/3V/ 2 I 1/4 („) 



+/3' sinhy 



+/3V/ 2 I 1/4 (y) 



4Jal(5v~ x 



\.^fa{(5v~ 



a\(5v 



-2a± v 



-3 



the following condition on the b p 

V b r(p + 1) f^Y /N = u vi\ 

p r(( P + i)/N)r((p + o/N) Km) 1 • ] 

Remarkably, this condition is already sufficient to cancel not only the leading exponential term but 
in fact the entire series of exponentially growing terms. Eliminating &at_2, the final solution for N 
integer becomes 

N-3 

Q(v) = bpfl p (v) 

p=0 

nf . PT t(p + 1) rgrji +jg) / ai \ (p+2-n)/n N _ 2 

ilJv) = v p s v n -j — — — - — ( ) v Fn-2 4.26 

py ' p r(£±i)r(2±i) r(iv-i) \n 2 J y 1 

Up to normalization, the form of Q(v) depends on £ and on N — 3 free parameters b p , while a\ merely 
sets a scale. The independent solutions fl p {p — 0, 1, . . . , N — 3) satisfy the boundary conditions 

Q p (v) ~ ( V " _ c ; V ~* ° (4.27) 



where explicitly 



/ cn\ CC+ri/iV T(^) r(p+l) 7rsin(f(p + 2)) 

^ r(i-c) r(^) r(^)sinte(p + c))sinfe(c-2)) [ ^ 

Therefore, we have not only eliminated the entire exponentially growing series, but furthermore, the 
fl p satisfy exactly the physically required boundary condition (see table |2|) for v — ► oo ]59 . Indeed, 
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Figure 3: Scaling functions v^Q p (v) for Typ I, N = 4, 



a i 



1 and p = 0,1. The different curves 



belong to values of £ = iVxi as indicated and the squares are the values of il } 



p,co • 



for N = 3 this cancellation of the exponential terms is a known property of the Kummer function 
iFi and we have for N = 3 the scaling function 



n 



N=3 



1^1 



— r 

°27T 



C 2 IT 



3 3 3ai 



r((c + i)/3)r(2/3) v 
r(c/3)r(4/3) (SaO 1 ^ 

C 2 t; 3 



i^i 



C 4 



3 3 3ai 



3 / V 3 3 3«i 



(4.29) 



where U is the Tricomi function and its known asymptotic behaviour U (a, b; z) ~ z~ a as 2 — > oo [|^, eq. 
(13.1.8)] reproduces the boundary condition eq. (|4.27| ). We have not found an analogous statement 
for N > 4 in the literature. Wright's formulas ||113|| simply list the dominant and the subdominant 
parts of the v —>■ oo asymptotic expansion but without any statement when the dominant part may 
cancel. Rather than giving a formal and lenghty proof of the cancellation of the entire asymptotic 
exponential series, we merely argue in favour of its plausibility through a few tests. In table ^ we 
list a few closed solutions for N = 4 with satisfy the boundary condition Q(0) = 1. By varying 
the parameters (3 and /3' one obtains the three independent solutions of the third-order differential 
equation ( 4.22|) . Furthermore, from the explicit form of the solutions we see that only the contribution 
parametrized by /3' diverges exponentially as v — > oo. There is a second solution which decays as 
for v — > oo and the third solution vanishes exponentially fast in the v — > oo limit. From the last 
two solutions we can therefore construct scaling functions with the physically expected asymptotic 



behaviour. In addition we illustrate the convergence of v^ p (v) towards Q P)00 , as given in (|4.28|) , by 



plotting v^Q p (v) as a function of v for several values of £■ This is done for iV = 4 in figure |3| and 
for N = 5 in figure ^ (we have also checked this for N — 6). Besides confirming the correctness of 
the asymptotic expressions fl4.27y4.28"D for v — » oo, we also see that for a large range of values of (, 
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Figure 4: Scaling functions v^Q p (v) for Typ I, A" 
( = Nxi. The circles indicate the values of f2 P)00 . 



5, q;i = 1 and p = 0,1, 2, for different values of 



the asymptotic regime is reached quite rapidly. 

Below, we shall need the explicit expressions for Q ^(v) for iV = 4 

r(3/4)^ r(n/2 + c/4) ( « 2 y 

n v ^ r((n+l + C)/4)s(n) 



^ = \/o F777TrT7ZY E 



2 r((C + l)/4) ^ r(n/4 + l)r((n + 3)/2) V </^i 

where s(n) := 4g (cos ^ + sin ^) cos ! | L . These expressions will be encountered again in section 5 
for the correlators of the ANNNI and ANNNS models at their Lifshitz points. 

Next, we study what happens for A" not an integer. It is useful to write the anisotropy exponent 

as 

l = N = N + e , e = V - (4.31) 
q 

where A" G N and p, q are positive coprime integers. 

For A" integer, we have seen that there is an unique solution which decays as Q(v) ~ t>~^ as 
v — > oo. The presence of such a solution for arbitrary A" may be checked by seeking solutions of the 
form 

oo 

Q(v) = ^a n v- n/q+s , a (4.32) 

n=0 

In making this ansatz, we concentrate on those solutions of eq. ( 4.22| ) which do not grow or vanish 
exponentially for v large. As done before for Typ II, and under the same conditions, we find upon 
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substitution 



oo^ aiT (l^ n+p + qN() ) + s+ l\ oo^, 1 . 

> 7 v -a n - P - q N v n/q + > -(n - px! - gA^oXi) - s a n t> " /9 = 



(4.33) 



which must be valid for all positive values of v. Comparing the coefficients of v, we obtain 

V 

s = — x x — N xi — — (N +e)xi — — ( 
Q 

a n = ; for n = 1,2, . . .,p + qN - 1 (4.34) 
n air(-i(n-p(l + xi)-gJV (l + xi)) + l) 

1 r(-^-^x 1 -N x 1 + 2) 

In principle, there might be additional ^-function terms which come from the definition ( |A2|) of the 
fractional derivative, see appendix A. However, if we either restrict to v > or else if ( is distinct 
from the discrete set of values ( c = 2 + m — n/q, where n, m G N, these terms do not occur. 

We can now let n = (p + qN )£ and a n = a( p+q N )e ='■ bi and find the simpler recurrence 

— oini+f-v+'M^ , <=1>2i ... (4 . 35) 



Ni r(2-(e + Xl )N) 

which can be solved in a manner analogous to the one used for Typ II before, with the result 

b =b ( ai Y r (i + (C-i)/AQ r(i-Q , , 

1 °\n 2 J e\r(i + i + ((-i)/N)r(i-(-£N) K ' } 

In the special case N = 1, the resulting series may be summed straightforwardly and leads to the 
elementary result &n=i{ v ) = b$v ^e _ai//l ', where v = tr~ 2 and ( = X\. We thus recover the form 
( p.24| ) of Schrodinger invariance for the two-point function (<f)<f)*), as it should be. 

If we use the identity r(x)T(l — x) = it/ sin(7rx) and define the function 

0<? W :=rmr(i- o) f Q^,< (437) 

the series solution with the requested behaviour Q(v) ~ boV~^ at v — *■ oo may be written as follows 



oo 

Q(v) = bj v 

£=0 

oo 



f v r (1 + (C - 1)/N) sin (ir(N£ + Q) T(£N + Q / a* J_ i 

^ r(C) sinvrC £!T(/ + l + (C-l)/iV) \N 2 v N I { J 



£=0 



From these expressions, it is clear that the radius of convergence of these series as a function of the 
variable 1/v is infinite for N < 2 and zero for N > 2. In the first case, we therefore have a convergent 
series for the scaling function, while in the second case, we have obtained an asymptotic expansion. 
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In several applications, notably the 3D ANNNI model to be discussed in the next section, the 
anisotropy exponent 9 ~ 1/2 to a very good approximation. Therefore, consider fractional derivatives 
of order a = Nq + e, where Nq is an integer and e is small. To study perturbatively the solutions for 
e -C 1, we use the identity ( |A12| ), set a = N = N + e and expand to first order in e. The result is 



rcg 



f iNo) (r) + eC f {No \r 

oo 

-(C E + lnr) + J2 



-1 



(4.39) 



t\t 



1 H 
- r 



where Ce = 0.5772 ... is Euler's constant. To find the first correction in e with respect to the solution 
( [4.26| ) when N is an integer, we set again N = N + e with iVo G N and e > and consider 

Q(v) = fi (0) (v) + en w (v) + O (e 2 ) (4.40) 

Then fi,(°\v) solves eq. ( 4.22|) with N = N and is consequently given by eq. ( [4.2(j| ), whereas Q^(v) 
satisfies the equation 

(a^ - 1 - v 2 d v - £v) n^(v) = oo(v) := -C {v 2 d v + (v) rt°\v) (4.41) 

which we now study. 

First, we consider the limiting behaviour of Q^(v) for v either very large or very small. If 



v > 1, we see from eq. (14.380 that f2 (0) ( 



iv) ~ v 



^(1 + 0(v N °))- This implies in turn that u(v) 



are some constants. Therefore one must have tt^'l 



(Aoo + B^lnv^-^ -^, where A^B^ 
v~**(l + 0(v~ N °)) in order to reproduce this result for u(v). On the other hand, if v -C 1, we 
have fi^(v) ~ cste. which leads to u(v) ~ (A + B lnv) with some constants A ,B . This can 
be reproduced from the limiting behaviour fiW(t>) ~ 0(f, f^ lnt> ). In conclusion, the first-order 
perturbation is compatible with the boundary condition ( |3.90p for the full scaling function Q(v). 

We now work out the first correction fiW(u) explicitly for N = 4. This is the case we shall need 
in section 5. There are two physically acceptable solutions &q\v) and fij (v) of zeroth order in e 
which are given in ( (4.30 ). The general zeroth-order solution is given by 

P 



nf ) («) + 



1/4 



(0), 



(4.42) 



where p = 6i«^ 4 /&o is a universal constant. Then all metric factors in the scaling function are 
absorbed into the argument v/ot^ 4 and the form of Q(v) is given by the two universal parameters ( 
and p. 

Consider the first-order correction to Q^°\v). From the explicit form of Cq and (|4.30 ), we have 



uo\v 



oo 

E 

n=0 



[A n v n+1 + B n v n+1 \nv] 



(4.43) 



where A n = —ip(n + l)B n and 



B,, 



T(n+C/4) r(3/4) 4n+C 
r(n+3/4) T(C/4) (2n)! 

1 y^r(n+(C+l)/4) 



1 

itti 



4n+l+C / 1 



Q i/4r((C+l)/4)r(2n+3/2) n!2 

_ 1 T(n+(C+2)/4) 4n+2+C ( J_ 
y/ai 2T(n+5/4) (2n+l)! ^4ai 





1^7 P 



n = mod 4 
1 mod 4 



i 



r(c/4) + Pr((c+i)/4) 



n 

n 
n 



(4.44) 



2 mod 4 

3 mod 4 
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Here ip(z) = T'(z)/T(z) is the digamma function [I]] and the identity 



a 



l+i 

— = C E + 4>{a + 1) (4.45) 



was used. The solution of the third-order differential equation ( 4.41 ) is of the form 



ftW ( v ) = J2 Kv n + b n v n In v] (4.46) 



71=0 



where in addition 
and, for all n > 



a = a 3 = , 6 = h = b 2 = b 3 = (4.47) 



ai [n + 4)(ra + 3j(n + 2) 

a n+4 = — ? — 1 — — [A n + a n (n + C) + &n- &n+4ai (3ri 2 + 18n + 26)1 (4.48) 
ai (n + 4)(n + 3)(n + 2) L v /J 

The value of the constant ao is fixed because of the boundary condition Q(0) = 1. The first pertur- 
bative correction fiW(u) still depends on the free parameters a\ =: p a and a 2 =: b. We also observe 
that because of ( 4.44J) , we have a^ n+3 = 64„+3 = for all n G N and that the metric factor ot\ merely 



v 



sets the scale in the variable v/a^ but does not otherwise affect the functional form of f^i 

We have seen above that for v ^> 1, we must recover Q(v) ~ v~^. We can therefore fix o and b 
such that the correction term fiW(w) goes to zero for v large. Furthermore, we see from figure^ that 
the asymptotic regime is already reached for quite small values of v for a large range of values of £. 
To a good approximation, we can therefore determine a and b from the requirement that QP'iyo) = 
and dQ^ (t> ) / dv — if v is finite, but chosen to be sufficiently large. The recursion (|4.48| ) then 
gives a system of two linear equations for a and b. 

As an example, we illustrate this in figure [5] for ( = 1.3. From figure |3], we observe that v o = 6 is 
already far in the asymptotic regime. For the values p = ±0.11, the scaling function (|4.40| ), with the 
first-order correction included, is shown for several values of e = N — 4. The first-order perturbative 
corrections with respect to the solution found for iV = 4 are quite substantial, even for small values 
of e. This suggests that a non-integer value of N in the differential equation ( |4.22| ) should be readily 
detectable in numerical simulations. We shall come back to this in section 5 in the context of the 
3D ANNNI model. 

The full series solution leads to difficulties with the boundary condition fi(0) = 1. This is further 
discussed in appendix D. 



5 Applications 

5.1 Uniaxial Lifshitz points 



Uniaxial Lifshitz points |)4] are paradigmatic examples of equilibrium spin systems with a strongly 
anisotropic critical behaviour. They are conveniently realized in spin systems with competing interac- 
tions. Besides the well-known uniaxially modulated magnets, alloys and ferroelectrics | |1 14| , |106| , [5TJ , 
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Figure 5: Perturbative scaling functions Q(v) eq. (|4.4Q|) for Typ I around N = 4, with ( = 1.3, 
«i = 1 and (a) p = +0.11 and (b) p = —0.11, for several values of e. 



recently found new examples include ferroelectric liquid crystals, uniaxial ferroelectrics, block copoly- 
mers, spin-Peierls and quantum systems ||109| , |111| , [7|, ^5], [L04j| . For the sake of notational simplicity, 



we merely consider uniaxial competing interactions, which are described by the Hamiltonian 



n 



SiS 



(5-1) 



where Si are the spin variables at site i. We shall consider here the ANNNI model,^ where Si = ±1 are 
Ising spins and the ANNNS model, where the Sj G M and satisfy the spherical constraint s\ = A/", 
where Af is the total number of sites of the lattice. The first sum runs only over pairs of nearest- 
neighbour sites of a hypercubic lattice in d = d± + 1 dimensions. In the second and third sums, 
additional interactions between second and third neighbours are added along a chosen axis (||) and 
en is the unit vector in this direction. Finally, J > and K\ t 2 are coupling constants. For reviews, 
see 



114, 106, 107, 81, 33 



In order to understand the physics of the model, we take k 2 = for a moment. If in addition K\ 
is small, the model undergoes at some T c = T c (k,i) a second-order phase transition which is in the 
Ising or spherical model universality class, respectively, for the systems considered here. However, 
if Ki is large and positive, the zero-temperature ground state may become spatially modulated and 
a rich phase diagram is obtained ||1 14|, |106|, |107|, jSl]. A particular multicritical point is the meeting 



point of the disordered paramagnetic, the ordered ferromagnetic and the ordered incommensurate 
phase. This point is called an uniaxial Lifshitz point (of first order) [Qj. If one now lets vary k 2 , one 
obtains a line of Lifshitz points of first order. This line terminates in a Lifshitz point of second order 
Lifshitz points of order L — l can be defined analogously and exist at non-zero temperatures 



83, 105 



for d > d* |[L05|| . For the ANNNS model, the lower critical dimension is 



4 := 2 + (L - I)/ L 



(5-2) 

6 The axial next-nearest neighbour Ising/spherical or ANNNI/S model is given by ( pTl| ) with K2 = 0. The case 
«2 7^ is sometimes referred to as the A3NNI model. For simplicity, we take here ANNNI/S to stand for axial 
non-nearest neighbour Ising/spherical and keep these abbreviations also for n-i ^ 0. 
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Close to a Lifshitz point, the scaling of the correlation functions C(r»,r±) is strongly direction- 
dependent. Here rn is the distance along the chosen axis with the competing interactions and r± is 
the distance vector in the remaining d± directions where only nearest-neighbour interactions exist. 
Slightly off criticality, correlations decay exponentially, but the scaling of the correlation lengths is 
direction-dependent 

en~(r-T L )-^ , £±~{T-T L )-^ (5.3) 

where Tl is the location of the Lifshitz point. The anisotropy between the axial (||) and the other 
(_L) directions is measured in terms of the anisotropy exponent 

= i/||/i/± (5.4) 

Precisely at the Lifshitz point, one expects 

C a (r h O)~r- 2x ° /e , a(0,r ± )~r^ , C £ (r||,0)~r^ /e , C £ (0, r ± ) ~ rj 2 - (5.5) 

for the connected spin-spin correlator C a and the connected energy-energy correlator C e , respectively, 
and where x a and x £ are scaling dimensions. The critical exponents a, (3, 7 are defined as usual from 
the specific heat, the order parameter and the susceptibility, but some of the familiar scaling relations 
valid for isotropic systems (where 6 = 1) must be replaced by 

2 - a = d ± u ± + 6v\\ , 7 = (2 - V± ) u ± = (2/6 - m ) n (5.6) 

where the anomalous dimensions m± are defined from the spin-spin correlator 

a(0,r ± )~r^ +e - 2+ ^ ] , CU^O)-^^ 2 ^ 1 ' 1 (5.7) 

and are related via m = rj±/9. Alternatively, one often works with exponents = v±, = u», 
V12 = i]± and rj u = m + 4- 2/6, see e.g. H§ |0|. Then 7 = (4 - r) ti )v i4 = (2 - % 2 )^2- 

Standard renormalization group arguments lead to the following anisotropic scaling of the corre- 
lation functions 

C„, e (r h r±) = b- 2x <^C ( ,, E {r ll b- e ,r ± b- 1 ) = r^n^ ( n /r e ± ) (5.8) 

for both the spin-spin and the energy-energy correlators, respectively. For a Lifshitz point in (d± + 1) 
dimensions, we have 

2(6 + d ± ) _ 2(g + d±)(l-a) r . 

W 6(2 + ^/(3) ' ^ £ 6(2 -a) { } 

We want to compare the form of the spin-spin correlator with the predictions of local scale 
invariance. We begin with Lifshitz points of first order. Then, as will be discussed further below, 
6 ~ I at least to a good approximation. In terms of the notation of sections 3 and 4, this corresponds 
to iV = 2/6 = 4. For iV = 4, we recall the two-point function of Typ I 

C(r ll ,r ± )=rJ_ <:/2 b (n ( v ) + J^Q^ , v = tr' 1 ' 2 (5.10) 

where VL ^(v) are explicitly given in eq. fl4.30|) . The functional form of Q(v) only depends on the 

universal parameters C and p. The metric factor ot\ only arises as a scale factor through the argument 

-1/4 

V OL x . 

We shall now present tests of the two-point function of Typ I of local scale invariance in three 
distinct universality classes. 
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1. Our first example is the exactly solvable ANNNS model. The phase diagram is well-known 
and uniaxial Lifshitz points of first order occur along the line P3] , |105| , [HJ 



«2 



1 , N 2 

1 — 4/ti) , K\ < 



9 v x/ ' A 5 
with a known Tj, = Tl(k\, Kq). The lower critical dimension 



(5.11) 

-. We need the following exactly 



known critical exponents in d dimensions [[83], |1 05 

4 



7 



2d -5 



2 



Co 



2U-- 



(5.12) 



which means iV = 4 in our notation. The exact spin-spin correlator along the line (|5.11|) of Lifshitz 
points is |40 



C*(r\\,r ± ) 



ty(a,x) 



d-d* 



II 



2 — 5ki rj_ 



fc=0 



T(k/2 + a) 



k\ T(k/2 + 3/4) 



x 



(5.13) 



(5.14) 



where Cq is a (known) normalization constant.^] This reproduces the exponent ( a from eq. ( |5.12 ). 
Comparing with the expected form (|5.10|) and the specific functions (|4.30|) , we see that 

,,2 



n (v) 



r ( 3 / 4 ) y {C 



r(C./4) 



4'2 



(5.15) 



With the correspondence t <-> ry, r «-> r± and the non-universal metric factor a% = |(2 — 5«i), we 
therefore observe complete agreement. In particular, we identify the universal parameter p = |j55| . 

2. Next, we consider the uniaxial Lifshitz point in the 3D ANNNI model. Two complimentary 
approaches have been used. First, the model may be formulated in terms of a n-component field 
(f>(r», r_|_) with a global 0(n)-symmetry and spatially anisotropic interactions [64]. This model, which 
might be called ANNNO(n) model, reduces to the ANNNI model in the special case n = 1 and gives 
the ANNNS model in the n 



oo limit. Recently, Diehl and Shpot 



3U 32, 33 studied 



8C 



very thoroughly the field-theoretic renormalization group of the ANNNO(n) model at the Lifshitz 
point at the two-loop level and derived the critical exponents to second order in the e-expansion, 
<i.0 Second, one may resort to numerical methods, such as series expansions | 



where e = 4.5 

or Monte Carlo simulations. While older simulational studies |68|] were restricted to small systems, 
the use of modern cluster algorithms [|112|[ allows to simulate considerably larger systems. The 
Wolff algorithm can be adapted to systems with competing interactions beyond nearest neighbours 
such as the ANNNI model [fXj |59|. In addition, a recently proposed scheme [Q permits the direct 
computation of two-point functions on an effectively infinite lattice. That technique can be extended 
to ANNNI models as well 



90, 59 



Before the scaling form of any correlator can be tested, the Lifshitz point must be located precisely. 
In table |4] we show some estimates for the coupling and the Lifshitz point critical temperature 
Tl- Here the ANNNI Hamiltonian (|5.1| ) with K2 = on a 3D simple cubic lattice was used. The 
increase in precision coming from the new cluster algorithm is evident and we take the estimates 
obtained in p0 as the location of the Lifshitz point. 



7 Properties of the function <3>(a,2;) are analysed in 108|. Explicit expressions are known for integer values of 
( = 4a and may be recovered as special cases of the functions listed in table [| 

8 Another recent two-loop calculation [^| apparently used some uncontrolled approximation in order to be able to 
evaluate the two-loop integrals analytically. See |3l], Q for a critical discussion. 
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Table 4: Estimates for the location of the uniaxial Lifshitz point of the 3D ANNNI model on a cubic 
lattice. The numbers in brackets give the uncertainty in the last digit (s). 



T L Kl,L 




3.73(3) 0.270(5) 
3.77(2) 0.265 
3.7475(50) 0.270(4) 


high-temperature series [36] 
Monte Carlo [Bfi] 
cluster Monte Carlo [|9Q|] 



Table 5: Estimates for critical exponents at the uniaxial Lifshitz point of the 3D ANNNI model. 
The numbers in brackets give the uncertainty in the last digit (s). 



a 


P 


7 


9 




0.20(15) 




1.62(12) 




high-temperature series 


8C 




0.19(2) 


1.40(6) 




Monte Carlo 


681 


0.160 


0.220 


1.399 


0.487 


renormalized field theory 


108, 


0.18(2) 


0.238(5) 


1.36(3) 




cluster Monte Carlo 


90 j 



Next, the anisotropy exponent 9 and the scaling dimension ( a must be found. While it had been 
believed for a long time that also for the ANNNI model 9 = ~ might hold, it has been recently 
established that to the second order in the e-expansion 9 = ~ — ae 2 + 0(e 3 ), where a ~ 0.0054 in the 
3D ANNNI model |HJ. In table |5| we list two older and the most recent estimates for the Lifshitz 
point critical exponents a,/3, 7 and 9. A direct determination of 9 from simulational data is not yet 
possible. Since in |Kj the exponents a, f3, 7 were determined independently, their agreement with 
the scaling relation a + 2/3 + 7 = 2 to within « 0.8% allows for an a posteriori check on the quality 
of the data. For details on the simulational methods we refer to fl90|, |59| . 



If we take the exponent estimates of |X| and in addition set 9 = 
ANNNI model (d ± = 2) 

Ca = 1.30 ±0.05 , Ce = 4.5 ±0.2 



2 ■ 



we find from (6.91) for the 3D 



(5.16) 



where the errors follow from the quoted uncertainties in the determination of the exponents a, (3, 7. 
If we now take 9 = 0.48 as suggested by two-loop results of ||108 , the resulting variation of both 
( a and ( £ stays within the error bars quoted in eq. ( |5.16| ). In conclusion, given the precision of the 
available exponent estimates, any effects of a possible deviation of 9 from | are not yet notable. We 
shall therefore undertake the subsequent analysis of the correlator by making the working hypothesis 
9=h 



In this case, we can compare with the scaling prediction (|5.10|) obtained for iV = 2/9 = 4. In 
figure 13 we show data PO] for the modified scaling function of the spin-spin correlator 



$(v) 



(5.17) 



The clear data collapse establishes scaling. It can be checked that there is no perceptible change in 
the scaling plot for values of 9 slightly less than | [90 |. 

For a quantitative comparison with Q5.10 ), one may consider the moments 



M(n) 



dv v n &(v) 



(5.18) 
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Figure 6: Scaling function as defined in ( 5.17|) for the 3D ANNNI model, for k,\ = 0.270 and 
T = 3.7475. Selected data on a 200 x 200 x 100 lattice are shown. The symbols correspond to several 
values of r±. The inset shows the full data set (gray points) and the prediction ( |5.10| ) of local scale 
invariance with N = 4 and p = —0.11, a± = 33.2 and b = 0.41. The data are from |0U . 



Table 6: Values of the parameters p and a\ as determined from different moment ratios 

Jk ({m} ;{nj})- 



k {™>i} { n j} 


P OL\ 


2 {0,-0.5} {-0.25,-0.25} 
2 {-0.25,-0.75} {-0.5,-0.5} 

2 {0.2,-0.9} {0,-0.7} 

3 {0.2,-0.6,-0.8} {-0.3,-0.4,-0.5} 
3 {-0.1,-0.6,-0.7} {-0.4,-0.5,-0.5} 


-0.102 32.7 
-0.125 34.0 
-0.100 32.8 
-0.102 32.8 
-0.117 33.5 



Then it is easy to show || that the moment ratios 

k I k k k 

J h ({m} ; {%}) = n M ( m *) / I1 M K) > with H m i = H n 3 ( 5 - 19 ) 

i=l I j=l i=l j=l 

and k > 2 are independent of bo and a\. They only depend on the functional form of $(w). This 
in turn is determined by ( a and p. Therefore, the determination of a certain moment ratio allows, 
with ( a given by ( |5.16D , to find a value for p. The Monte Carlo data will be consistent with (|5.10| ) if 



the values of p found from several distinct ratios coincide. In practise, these integrals cannot be 
calculated up to v = oo but only to some finite value vq and the moments retain a dependence on 
ct\ through the upper limit of integration. Then an iteration procedure must be used to find p and 
a± simultaneously [pOj] . The results are collected in table |6|. 



Clearly, the two parameters can be consistently determined from different moment ratios. The 
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final estimate is [|90 



p = -0.11 ±0.01 , ai = 33.2 ±0.8 (5.20) 

Since we have seen above that p = for the spin-spin correlator of the ANNNS model, it follows 
that the value of p is characteristic for the universality class at hand. 

In figure ^| the Monte Carlo data are compared with the resulting scaling function, after fixing 
the overall normalization constant b = 0.41. The agreement between the data and the prediction 
( p.lUj ) of local scale invariance is remarkable. 

To finish, we reconsider our working hypothesis 9 — \. Indeed, in figure |^ we had shown how 
the form of the scaling function Q(v ) changes when e = N — 4 is increased, to first order in e. In 
particular, rather pronounced non-monotonic behaviour is seen for values of e ~ 0.1 which is the 

see table ||. 



order of magnitude suggested from the results of renormalized field theory |3Q, 108 



Nothing of this is visible in the Monte Carlo data of figure || Assuming that first-order perturbation 
theory in e as described in section 4 is applicable here, we conclude from this observation that e 
should be significantly smaller. Given the differences between the exponent estimates coming from 
renormalized field theory [ |108| | and cluster Monte Carlo [90|, a possible difference of 9 from | cannot 
yet be unambigously detected. Direct precise estimates of 9 are needed. 

A similar analysis can be performed for the energy-energy correlation function. This will be 
described elsewhere. 

In summary, having confirmed local scale invariance for the spin-spin correlator at the Lifshitz 
points in the ANNNI and the ANNNS models, it is plausible that the same will hold true for all 
ANNNO(n) models with 1 < n < oo. 

3. Finally, we consider the Lifshitz point of second order in the ANNNS model. In the ANNNS 
model as defined in eq. (|5.1| ), a second-order Lifshitz point occurs at the endpoint 

2 1 

5 ' K2 ~ ~15 



(5.21) 



of the line (|5 . 1 1|) [|105|1 . The lower critical dimension = |. We need the following critical exponents 



105] 



p= 1 - 

H 2 



7 



6 



3d -8 



8 = 1 
3 



Co 



d-- 

3 



(5.22) 



which in our notation corresponds to iV = 6. Therefore, the prediction of local scale invariance is 
G(t,r) = r~^ //3 fi(tr~ 1//3 ) with £l(v) given by Q4.23j , f05|) . At the Lifshitz point, the exact spin-spin 
correlation function is |4(J (with c = 144/5) 



1 



1 



-1/3 I 1/3 



(5.23) 



where the scaling function^ can be written^ in terms of generalized hypergeometric functions 1-F4 



Via] 



7rr(5/6) 
3T(a± 1/3) 



1F4 



112 5 
l ' 3'2'3'6'" 
1 2 5 7 4 



x 

27 



XlFi V a+ 3 ; 3'6'6'3 ; 27 J H ^T(l/6) 



6T(a±2/3) 



x 2 1-F4 ( a + 



(5.24) 
2 7 4 3 5 x 3 \ 



9 Properties of the function 5(3, a; x) are analysed in f 
for a = n + 5, n + I with n 6 N. 

10 We correct herewith a typographical error in eq. (4.1) in Eol 



3 6 3 2 3 27/ 
Explicit expressions in terms of Airy functions are known 
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and Co is a normalization constant. This reproduces the exponent ( a from eq. (|5.22 ). We now 
compare the function S(3, a; x) with the expected form ( |4.23[ ) with N = 6. The arguments of the 
functions S and fl are related via x 3 = — A— which implies 

48 ot\ 1 



,i7T/3 



a; 



(48ai) 1 /3 



(5.25) 



Using this correspondence, we get 



c 



— I 3, ; x 



r(a) 



vrr(5/6) 



3T(a + l/3) e 



m/3 



(48a 1 ) 1 /3 



2 6T(a + 2/3) e 

V J- 2 + 



2i7r/3 



iT(l/6) (48ai) 2 /3 



(5.26) 



where the T v are defined in eq. ( |4.23| ). The general form of the scaling function for = 6 is 
= Yl A p=QbpV p J r p and we can identify the values of the free parameters b v which apply for the 
spin-spin correlator at the Lifshitz point of second order in the ANNNS model. We find b\ = 63 = 
and 



r(o) 



3T(a + l/3) e 



m/3 



6T(a + 2/3) e 



2m/ 3 



(5.27) 



/iT(5/6) ' "* 7T (48« 1 ) 1 /3 ' 4 v^r(l/6) (48a!) 2 / 3 

It is now straightforward to check that the constraint eq. ( |4.25| ) is indeed satisfied. In view of the 

1 ■ 

2' 



known [40 power-law decay of the function S(3, a;x) for x — > 00 (and a ^ |, |) this result might 



have been anticipated. 

We are not aware of any study of a second-order Lifshitz point in a different model. 



5.2 Aging in simple spin systems 

We now turn to a class of systems which display non-equilibrium dynamical scaling. For the sake 
of simplicity, consider a simple ferromagnetic spin system, e.g. an Ising model, evolving according 
to some dynamical rule. We shall exclusively consider the case of a non-conserved order parameter. 
Prepare the system in some initial state (an infinite-temperature initial state without any correlations 
is common) and then quench it to some fixed temperature T below or equal to the equilibrium critical 
temperature T c . Then follow the evolution of the system at that fixed temperature T. In the first case, 
the system undergoes phase-ordering kinetics while in the second case one considers n on- equilibrium 
critical dynamics. In both cases, the equilibrium state is never reached for the spatially infinite 
system. Rather, correlated domains of typical time- dependent size L(t) ~ t 1 '" form and grow, where 
z is the dynamical exponent. Consequently, the slow dynamics displays several characteristic features 



which are absent from systems in thermodynamic equilibrium. For reviews, see |L3], y| [20|, |48| . 

Main observables are the two-time correlation function C(t, s; r) and the two-time response func- 
tion R(t, s; r), defined as 



C{t, S] r-r') = (a r (t)a r ,{s)) , R{t,s;r-r r ) (){a ' A ' h) 



5h r >(s) 



(5.28) 



where o~ r (t) is an (Ising) spin variable and h r (t) the conjugate magnetic field at time t and at the site r. 
It is assumed throughout that the quench occurred at time zero and that spatial translation invariance 
holds. In particular, we shall focus here on the two-time autocorrelation function C(t, s) = C(t, s; 0) 
and the two-time autoresponse function R(t,s) = R(t,s;0). Then s is the waiting time and t the 
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observation time. If either T < T c or T = T c the system is always out of equilibrium in the sense 
that the fluctuation-dissipation ratio [f2^, [27], [T5| 



X{t,s) = TR(t,s] 



( dC(t, 
V ds 



7^1 



(5.29) 



Furthermore the two-time observables such as C = C(t,s) and R = R(t,s) depend on both the 
waiting time s and the observation time t and not merely on their difference r = t — s. This breaking 
of time translation invariance is usally referred to as aging [15], |20], |48| and will be used in this 
sense from now on. In addition it is well-established |13| that the aging process is associated with 
dynamical scaling, that is in the scaling limit s — > oo and t — > oo such that 



x 



t/s > 1 



is kept fixed, one has 



C(t,s) 



■,-b 



fc(t/s) , R(t,s) 



-l-a 



fn(t/s) 



(5.30) 



(5.31) 



where a, b are non-equilibrium critical exponents and fc and are scaling functions. For large 
arguments x ^> 1, these scaling functions typically behave as 



fc(x) 



x 



-X c /z 



f R {x) 



X 



(5.32) 



where Ac, \r are the autocorrelation |38|, |65|] and autoresponse exponents.^ Remarkably, it can be 
shown that at late times the form of the growth law L = Lit) (and thus the value of the dynamical 
exponent z) can be found for phase-ordering kinetics of purely dissipative systems from the scaling 
of the two-time correlation function C(t, s; r) p5| . 

A. If in 



one has Ac = A 



For fully disordered initial conditions, recently reviewed in [4S 
addition one has T = T c , the relation a = b = 2(3 / vz holds below the upper critical dimension, where 
(3, v are standard equilibrium critical exponents and the critical autocorrelation exponent A = d — z<3, 
where G is the initial-slip critical exponent f66| . If on the other hand T < T c , one has 6 = 0, but 
there does not seem to exist a general result for a. Indeed, in the Glauber-Ising model a = 1/2 in 
2D and in 3D, while in the kinetic spherical model a = d/2 — 1. However, these exponent identities 



do not necessarily hold for more general initial conditions [0, pq| . We shall need here the values of 
the exponents z, a and Xr which are collected in tables [7| and || below, for the Glauber-Ising model 
and the kinetic spherical model, respectively. 

We now derive the exact form of the scaling function fn(x) for the autoresponse function and 
then generalize towards the full spatio-temporal response function R(t,s;r). Afterwards, we shall 
describe tests of these predictions in specific models. 

We begin by assuming that the response functions transform covariantly under local scale trans- 
formations flIB| . Recall that in the context of Martin-Siggia-Rose theory, see 
therein, these are given in terms of correlators 



191 and references 



R<M,{ti,t 2 ]r 1 ,r 2 ) 



6hW(t 2 ,r 2 ) 



(0(ti,ri)^(t 2 ,r 2 )> 



(5.33) 



/i«0=0 



of the scaling operator <p{t, r) and the response operator ip(s, r) associated with the field h^' canon- 
ically conjugate to the scaling operator ip. Usually, one merely considers the response function 

11 The values of the exponents Ac, (and also a, b, z) depend on whether T < T c or T = T c , but we shall use the 
same notation in both cases. 
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R = R^ = (00) eq. ( |5.28| ) of the order parameter with respect to its own conjugate magnetic field. 



There is a clear analogy to the correlators (00*) which we have found in sections 3 and 4. 

However, the treatment carried out in sections 3 and 4 cannot be entirely taken over since in aging 
systems time translation invariance is broken. Therefore, the autoresponse functions R = R(t, s) does 
not transform covariantly under the entire set of local scale transformations constructed in sections 



3 and 4 but only under those belonging to the subalgebra [53 



S = {X ,X 1 ,Y m ,...} (5.34) 

This contains both scale (Xq) and special conformal (X\) transformations as well as space translations 
(YLiAj) together with all generators which are obtained from the commutators of these, see section 3. 
Since the dynamics of aging systems may thought of as being described by some Langevin equation 
which is of first order in time, the realization of Typ II of local scale transformations will be adequate. 
From the explicit form of the generators eq. ( 3.72j) we see that the line t = is kept invariant under 



the action of the X 01 and the Y m . 

We therefore consider a two-point function R = R(ti, fi, ^2) of two scaling operators 0i t 2 
which transform covariantly under the action of S. These are characterized in terms of their scaling 
dimensions X{ and the parameters ft,7i, with % = 1,2. The covariance of R is expressed by the 
conditions ( |3.83 ) restricted to those generators contained in S. Now, spatial translation invariance 



can be implemented as shown in section 3 and leads to R = R(ti,t 2 ; r) with r = r% — r 2 provided 
only that the constraints 

ft + (-l) 2 ~ z Pi = , 72 + (-l) 2 -^7i = (5.35) 

hold, in complete analogy with eq. ( |3.91| ). 

The autoreponse function R = R(t, s) is now obtained by setting r = \r\ = 0. Then the conditions 
Y m R = are automatically satisfied and the last two remaining conditions X R = X\R = lead to 
the following differential equations 



(td t + sd s + — + — ) R(t, s) = 

\ z z J 



t 2 d t + s 2 d s + —t + —^Rii.s) = () (730) 



with the solution 



R(t, s) = r [ 2 ° (t - s)- {xi+X2)/z 0(t - s) (5.37) 



where r$ = ro(/3i, ft, 7i> 72) is a normalization constant which vanishes if the constraint ( 5.35 ) is 
not satisfied. We have also explicitly included the function which is required because of causality 
(Q(x) — 1 if x > and Q(x) =0 otherwise), see e.g. fn], [19|| . We can now compare with the 



expected scaling form eq. ( |5.31| , |5T32| ) and then arrive at the final result 

(, \ l+a— Xr/z 
- (t-s)- 1 ^ e(t-s) (5.38) 

i.e. f R (x) = r x 1+a - XR/z {x- ly^ Qix- 1) 

Therefore, once the exponents a and Xr/z are known, the functional form of the two-time autore- 
sponse function is completely determined. 
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We note the following interesting consequence: from the constraint ( |5.35| ), we see that for the 
response of the scaling operator to a perturbation by its own conjugate field to be non-vanishing 
one must have l/^l = |/3t| and \ = For different scaling operators (f> ^ i/j, the absolute values 
of ^ and (3^ will a priori be different (and similarly for | 7( ^| and |7^|) and therefore we generically 
expect Rfaj) = 0. 

Next, we derive the spatio-temporal response function R = R(ti,t 2 ;r) from its covariance under 
the action of S. It is convenient to write R in the form 

/+ \ ( x 2~Xl)/z 

R= (j-) G ( t ' r ) ' * = *i -*a (5.39) 

where we already took the explicit solution ( |5.37| ) into account. Since the generators Y m (with 
m = — 1/z, — 1/z + 1, . . .) do not modify the temporal variables, the treatment of the conditions 
Y m G = of section 3 goes through and we merely have to consider the two covariance conditions 
X R = XiR = 0. These lead to the differential equations for G = G(ti — t 2 ,r) (written here for 
d = 1 spatial dimensions) 

-td t — -rd r (xi + x 2 ) ) G = 

z z J 

. t % _ l trdr _ - z ^^t - (fix + 7l ) r 2 d 2 r ~ z - 7 i2(2 - zyd 1 -^ G (5.40) 

+2t 2 X G + -r 2 Y„ 1/z+1 G = 

z 

Therefore, the function G = G(t, r) does satisfy exactly the same equations (|3.92|) which were obtained 
in section 3 for a two-point function of quasiprimary operators 0^2 with effective scaling dimensions 

%l,e& = x 2 , c s = - {xi + x 2 ) (5.41) 

The solution of these equations may therefore be taken from section 4 and our final result for the 
space-time reponse is 

R(t, s; r) = R(t, s) $ ( _| r J ) (5.42) 

where the autoresponse function R(t,s) is given in eq. ( |5.38| ) and the scaling function $(n) can be 
read from eq. ( |4.16|) 

[ E zA (-2z(2 - z)7iu z ) ; 7i = -A 

where A = (z — l)/z + 2(2 — z)/[z(l + /?i/ 7 i)] is a universal constant and the functions (£ a ,&(w) and 
E a>b {u) are defined in eq. ( |4 . 1 8| ) . Eqs. ( |5.38| , |5.4^ ) are the main results of this section. 

For the special case z = 2 eq. (|5.42|) takes the simple form |53| (see also eq. ( |2.26|) ) 

/ M r 2 \ 

R(t, s; r) = R(t, s) exp (- — —] (5.44) 

where M. = 2(/?i + 7 i) is a dimensionful and non-universal scale factor. 

For simplicity of notation, the derivation has been carried out in the case of one spatial dimension, 
d — 1. If spatial rotation invariance holds, our result ( (5.42j ) also holds for d > 1. Indeed, as they 
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stand, rotation invariance is tacitly assumed in eqs. (|5.42 ) and (|5.44| ). However, if rotation invariance 



in space is broken (and in fact this is argued to be the case in phase-ordering kinetics [^B], 0), this 
means in our framework that the above calculation must be carried out separately in each spatial 
direction. Consequently, eq. ( p.42[ ) still holds phenomenologically, but the non-universal constants 
/?i,7i should become direction-dependent. 

We now turn towards tests of the conformal invariance prediction (|5.38| ) in specific models. Either 
experimentally or in simulations, it is hard to measure R(t, s) directly and one rather studies the 
integrated response. For definiteness, we shall use from now on the Ising model language, since 
afterwards the numerical tests will be performed in that model. Depending on the history one 
obtains different forms of the integrated response. For example one may obtain the thermoremanent 
magnetization Mtrm(£, s) by quenching the system in a small magnetic field h, kept constant between 
the quench at time zero and the waiting time s and subsequently switched off. Alternatively, one 
may quench in zero magnetic field, switch it on at the waiting time s and keep it until the observation 
time t when the zero-field cooled magnetization MzFc{t, s) is measured. We then have two integrated 
reponse functions 



p(t,s) = T [ duR(t,u) = ^M TRM (t,s) 
Jo h 



'o 

X{t,s) = T ^duR(t,u) = ^M ZFC (t,s) (5.45) 

(often, p{t, s) is also called a relaxation function and x(^, s) a susceptibility function). Other histories 
are possible. In practice, for ferromagnetic systems h is a random magnetic field with zero mean in 
order to treat all phases equally |J. The scaling of the TRM integrated response is readily found 
from eq. ( |5.37|) 



p(t, s) = r T s i-(«-h»)/* x -^i/m Fl ( ! + ^_J^ 2 + Tl^l. X ~A (5.46) 

\ z z z J 

where 2F1 is a hypergeometric function and x = t/s. In terms of the exponents a and Xr/z this 
becomes 

p{t, s) = r T s~ a x~ Xr/z 2 F 1 (l + a, ^ - a; ^ - a + 1; x~^ (5.47) 

Once the values of the exponents a and Xr/z are known, the functional form of p(t, s) is completely 
fixed (the typical behaviour of the scaling function s a p(t, s) is exemplified in figures and [] below). 
Evidently, x(t, s) = pit, t) - p(t, s). 

We are now ready to compare these predictions with simulational and exact results in specific 
models. 



1. First we consider the kinetic Ising model with Glauber dynamics in both 2D and 3D p5| . The 
spin Hamiltonian is Ti = — Yluj) <T i <T j with Ising spins o~i = ±1 at site i and where the sum is over 



the nearest-neighbour pairs of a hypercubic lattice. Glauber or heat-bath dynamics is realized 
through the stochastic rule Si(t) — > Si(t + 1) such that 

8i(t + 1) = ±1 with probability |[1 ± tanh(/i i (t)/T)]. (5.48) 

The local field acting on Sj is hi(t) = h + 2^j(i) s i(^)) where h is the external magnetic field, 
denotes the nearest neighbours of the site i. The TRM integrated response p(t, s) was measured |p8| , 
following the method of p], B?J and using a small random magnetic field, for systems with 300 x 300 



46 



Table 7: Critical temperature and some non-equilibrium exponents of the 2D and 3D Glauber-Ising 
model with an infinite-temperature initial state, both for the phase-odrering regime (T = 0) and for 
critical dynamics (T = T c ). In 2D, T c = 2/ ln(l + yg) exactly. 





2D 


3D 


T 

A c 


2.2692 


4.5115 


Z 


T = 
T = T C 


2 

2.17 


2 

2.04 


Xr 


T = 

T = T C 


1.25 
1.59 


1.50 
2.78 


a 


T = 
T = T C 


0.5 
0.115 


0.5 

0.5064 



spins in 2D and 50 x 50 x 50 spins in 3D, with a fully disordered initial state. Larger systems 
were also simulated, in order to check for finite-size effects and averages over at least 1000 different 
realizations of the systems were performed. For a comparison with eq. (|5.47|) , which is equivalent to 
a test of (|5.38p , the relevant exponents are collected in table [7|, see |[8], ^5], p3 



Clearly, from ( |5.47| ) we expect a data collapse if s a p{t, s) is plotted against x = t/s. In figure [7| 
Monte Carlo data at criticality are shown and scaling is indeed seen to hold. Furthermore, upon 
adjusting the normalization ro, there is complete agreement between the data and eq. (|5.47|) . Sim- 
ilarly, data for T < T c are shown in figure [8|. While the expected scaling of s 1//2 p(t, s) as a function 
of x = t/s works well in 3D, that is not the case in two dimensions. However, recall that analyti- 
cal calculations [|10j on the scaling of the response function in the spirit of the OJK approximation 
rather suggest in two dimensions the presence of logarithmic corrections p(t,s) ~ s~ 1//2 In (s) /(t/s). 
Therefore, the following ansatz for the 2D Glauber-Ising model [pH[ 



p(t,s) 



,-V2 



ro + r\ In s) E ( - 



(5.49) 



appears natural, where r^i are non-universal constants and E = E(x) is a scaling function.0 Indeed, 
a satisfactory scaling is found this way, as can be seen from figure [|. It has been checked that the 



same kind of scaling holds in the entire low-temperature phase, where ro and r\ depend on T [p8 



Furthermore, the form of the scaling functions thus obtained are again in perfect agreement with the 
prediction eq. ( |5.47| ) of local scale invariance, for both 2D and 3-D. 

2. Second, eq. ( |5.38| ) has been tested extensively in the exactly solvable kinetic spherical model 
with a non-conserved order parameter, for d > 2 space dimensions. The kinetic spherical model may 

81 ffn PT7L [171, E3L 



be introduced as a spin model 
the J, 



with Hamiltonian 7i 



JijSiSj, where 



h0 are coupling constants and the are real variables subject to the spherical constraint 



£^ 2 = a/- 



(5.50) 



where the sum runs over the entire lattice and M is the number of sites. The dynamics of the model 
is generated through a stochastic Langevin equation 

dS r 5H[S\ 



dt 



- (2d + $(t))S r + 7lr(t) 



(5.51) 



12 It is not impossible that the logarithm might be explained through logarithmic Schrodinger invariance, which 
would have to be constructed by analogy with logarithmic conformal field theories, see e.g. ]39|, We hope to come 
back to this in the future. 
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Figure 7: Scaling of the TRM integrated response function p for the 2D (above) and the 3D (below) 
Glauber-Ising model at criticality (T = T c ) . The symbols correspond to different waiting times. The 
full curve is the local scale invariance prediction Q5.47]) for p(t, s). The data are from 



where the Gaussian white noise rj r (t) has the correlation 

(Vr(t)Vv'{t')) = 2T6 r yS(t - t') 



(5.52) 



and 3(t) is determined by satisfying the spherical constraint. Alternatively, the same model may be 
studied as the n — * oo limit of the coarse-grained O(n) vector model using field theory methods, see 



12, 23, pi 82, 16 



In many studies, the short-ranged spherical model with only nearest-neighbour interactions Jj j 
and an infinite-temperature initial state without any correlations was considered 52, [28|, ^7|, |117| , 
24[ . Exact results for the two-time autocorrelators and autoresponse functions were found. Writing 
the two-time autoreponse function in the scaling limit as R(t, s) = (An s)~ d l 2 f u(t / s) , one has in the 
ordered phase (T < T c ) |82 , 



(here and below, we always take t > s or x > V 
f R (x) = (x - l)- d / 2 



for all values of d > 2. At the critical point (T = T c ) one has [66, [47 1 

l)-d/2 



x l-d/if x 



x 



iyd/2 



if 2 < d < 4 
if 4 < d 



(5.53) 



(5.54) 
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Figure 8: Scaling of the TRM integrated response function p for the low-temperature Glauber-Ising 
model in 2D (above) at T = 1.5 and in 3D at T = 3 (below). The symbols correspond to different 
waiting times. The full curve is obtained from (|5.47|) . The data are from ]58 . 



and both results are in full agreement, upon identification of exponents, with eq. ( |5.38| ). The second 
expression of ( |5.54|) corresponds to the mean-field case and coincides with the result found for a free 
Gaussian field as expected. We point out that the form of the correlation scaling functions 
depends strongly on d and on whether T < T c or T = T c since the behaviour of the fluctuation- 



dissipation ratio X(t, s) is different [4~7] in each of the three cases considered so far. 

More recently, these results have been generalized in two directions. First, the two-time auto- 
correlation and autoresponse functions were calculated exactly in the kinetic spherical model with 
spatially long-range interactions, of the form [JF7] 



J, 



hi 



J{r 



hi / 



Kl I 



(5.55) 



where ^ ■ runs over all sites j ^ i, is the distance between sites i and j, Jo is a constant and a is a 
free parameter. For d > 2 and o > 2, one recovers the short-ranged spherical model discussed above. 
On the other hand, if either (i) d > 2 and < a < 2 or else (ii) d < 2 and < a < d the model has 
an equilibrium phase transition at a non- vanishing T c between an ordered and a paramagnetic phase 
In this case and below criticality (T < T c ), the dynamical exponent z = a and the response 
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Table 8: Some non-equilibrium exponents of the kinetic spherical model with short-ranged inter- 
actions and correlated initial conditions of the form ( p. 57] ), for the five distinct regimes I,. . . ,V at 
criticality (T = T c ) and also in the ordered phase (T < T c ). Here D = d + a + 2. 



Regime 


conditions 


F 


a 




z 


\r 


I 


2 < d < A 


2 < D < 4 


a/4-1/2 


d/2- 


1 


2 


d- a/2-1 


II 


4 < d 


2 < D < 4 


(d + a)/4- 1/2 


d/2- 


1 


2 


(d-a)/2 + l 


III 


2 < d < 4 


4 < D 


l-d/4 


d/2- 


1 


2 


3d/2-2 


IV 


4 < d 


4 < D a > -2 





d/2- 


1 


2 


d 


V 


4 < d 


4 < D a < -2 





d/2- 


1 


2 


d 


T <T C 


2<d 


(d+a)/4 


d/2- 


1 


2 


(d-a)/2 



function is in the aging scaling limit ]1 

R(t,s 



r (t/ s 



,d/(2a) 



(t - s 



-d/a 



(5.56) 



which again fully confirms eq. ( |5.38| ). One identifies the exponents a = d/a — 1 and \r = d/2 in the 
ordered phase. We are not aware of any published results on R in this model at criticality. 

Second, the case of nearest-neighbour interactions but with long-ranged correlations characterized 
by the form 

CiBi(T-) ~ \r\~ d ~ a (5.57) 

of the spin-spin correlator in the initial state has been studied [|S2"| , |T!| [23| , |88| . The uncorrelated 
initial state considered above is recovered as the special case a = 0. Indeed, the effect of initial 
correlations is only notable in the long-time behaviour if a < 0. In the ordered phase (T < T c ), the 
exact autoresponse function is in the scaling limit R(t, s) = (Aix s)~ d l 2 f 'jt(t / ' s) , where [^, |88|1 



x (d+a)/4 f 



(5.58 



in complete agreement with ( |5.38| ). For a = 0, eq. (|5.53|) is reproduced. At criticality (T = T c ), five 
distinct regimes of non-equilibrium critical dynamics exist p8| . These are distinguished by the values 
of d and a as listed in table § which also gives the values of the required non-equilibrium exponents. 

-d/2. 



The autoresponse function is R(t, s) 



- (An)' 
fii(x) 



1 a /j?(^/ s ) i n au fi ye regimes where 

l)-d/2 



x [x 



(5.59) 



and the values of the exponent F can be read off from table |^. Once more we find complete agreement 
with eq. ( |5.38|) . The results quoted in eq. (|5.54 ) are reproduced in the critical regimes III and IV. 
Again, the behaviour of the autocorrelations differs widely between the five regimes as discussed in 
detail in 



3. Recently, the off-equilibrium response and correlation functions of the 0(n)-symmetric vec- 
tor model with a non-conserved order parameter (model A in the terminology of [51]) and a fully 
disordered initial state were calculated in one- loop order in d = 4 — e dimensions [pl| . From earlier 
calculations of correlators and responses with the initial state at criticality |66], the critical exponents 



0(, 2 ) 



d 



-1 + {e 2 ) 



A 



R 



. e n + 2 , 2N 
d - + O (e z ) 



2 K ' 1 1 2n+8 

were already known. At criticality, the autoresponse function is in the scaling limit JIB 



R(t, 



r 



€ n + 2 
4 n+8 



(t - s y d/2 + O (e 2 ) 



(5.60) 



(5.61) 
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which again agrees, to first order in e, with the prediction ( |5.38| ). It would be extremely interesting 
to see whether this agreement can be extended to higher orders in the e-expansion. We are not aware 
of any results on R(t, s) in the low-temperature phase. 

4. A particularly simple textbook system which reproduces eq. ( |5.38|) is the free random walk, 
2D]. It is described in the continuum through the following Langevin equation for y(t) and the 



sec 



correlator of the Gaussian white noise r\ 

dy(t) 



dt 



r,(t) , ( V (t) v (s))=2T8(t~s) 



(5.62) 



The autocorrelation function and the autoresponse function are readily found [2c 



C(t,s) 
R(t, s) 



(y(t)y(s)) 
S(y(t)) 



5h(s) 



h=0 



2Tmin(t, s) 



(5.63) 



and do satisfy the scaling fom ( |5.32j ). The system is out of equilibrium since the fluctuation- 
dissipation ratio X(t, s) = 1/2. The exact expression for R(t, s) matches eq. ( |5.38j ) with the exponents 
a = — 1 and Xr/z = 0. 



5. Having studied so far simple ferromagnets, we now consider the spherical model spin glass [28 
as a very simple example of a disordered system. The Hamiltonian H = — Yluj) ^hj^i^j describes 
nearest-neighbour interactions between the spherical spins subject to the constraint ( [5.501 ). The 
couplings Jij are independently distributed quenched random variables with zero mean and variance 
inversely proportional to the number of sites N ' . For uniform initial conditions, the mean-field two- 
time autoresponse function is in the scaling limit [EH] 



R(t, s) = r 



t\ 3/ \ X 

s) 



■3/2 



(5.64) 



which again agrees with ( |5.38| ). Could this be a hint that the form ( |5.38| ) for R(t,s) might also 
hold for glassy systems ? In fact, the result ( |5.64j) coincides with eq. ( |5.54j) for d = 3. This is a 
consequence of the known ||117|| similarity between the 3D spherical ferromagnet and the mean-field 



spherical spin glass. Therefore, although the result ( 5.64 ) may appear suggestive, it is not yet clear 
at all whether or not the predictions ( |5.38| , |5"~4"7| ) of local scale invariance may be reproducible in 
physically interesting glassy systems. 

6. The tests described so far all concerned the autoresponse function R(t,s). Spatio-temporal 
responses R(t, s; r) have so far only been tested in the exactly solvable spherical model. For short- 
ranged interactions, the dynamical exponent z = 2 and from local scale invariance (or Schrodinger 
invariance in this case) we expect ( |5.44[ ) to hold. Indeed, this had been checked long ago for a 



disordered initial state both at and below T c |66|, |82], Q and recently for an initial state with the 

J. To leading order in the 



long-range correlations of the form ( |5.57|) , again both at and below T c 
^-expansion, eq. ( |5.44j ) also holds in the 0(n)-model |IBJ. Tests of the spatio-temporal response 
( |5.42j ) in the Glauber-Ising model both below and at criticality are currently being performed and 
will be reported elsewhere. 



Summarising, we have seen that the form eq. ( |5.38| ) of the two-time autoreponse function has 
been reproduced in a considerable variety of non-equilibrium ferromagnetic spin systems with a 
non-conserved order parameter. The examples for which ( |5.38[ ) has so far be seen to hold suggest 
its validity independently of the following characteristics of specific models, namely (see also note 
added) 
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1. the value of the dynamical exponent z. 

2. the value of the space dimension, provided d > 1. 

3. the number of components of the order parameter and the global symmetry group. 

4. the spatial range of the spin interactions. 

5. the presence of spatial long-range correlations in the initial state. 

6. the value of T, in particular whether T < T c or T = T c , provided T c > 0. 

The response functions of different universality classes are only distinguished by the values of the 
exponents a and Xr/z. 

Clearly, further tests in different models will be most useful to either confirm further or else 
invalidate this conjecture. On the other hand, as yet there are only few tests of the full spatio- 
temporal response (|5.42 ). At present, there is no prediction from local scale invariance for the 



two-time correlation functions in non-equilibrium dynamical scaling. 

Finally, we comment on some examples where eq. ( |5.38|) does not hold. From the exact solution of 
the ID Glauber-Ising model at T = one has in the aging regime R(t, s) = [2Tc 2 s(t — s)] -1 / 2 [0, 76| . 



For the 2D XY model in the critical low-temperature phase, one has in the spin-wave approximation 
R(t,s) ~ s^-^fnd/s) with f R (x) = [(x + lf/x\^ A {x - l)- 1 "'/ 2 0. Here rj = rj(T) is the 



usual equilibrium temperature-dependent exponent, see |0, |34], |116| , [[9] and references therein. It 



is well-established, however, that in these models already the growth law L = L(t) is unusual: the 
generic description which assumes that energy dissipation is dominated by the motion of single defect 
structures at the scale L = L(t) no longer applies here fl95 , 14| . From that perspective, it is not too 



surprising that also the form of the scaling functions should be non-generic in these models. We shall 
come back elsewhere to the question whether these systems may or may not be treated by some form 
of local scale invariance . 



5.3 Equilibrium critical dynamics 

Having discussed the local scaling of the response function out of equilibrium, we now briefly turn 
towards the case of equilibrium critical dynamics. It will be of interest to compare the predictions of 
local scale invariance as introduced here with those obtained from dynamical conformal invariance 



1< 



Consider a spin system at its equilibrium critical point. Two-time correlation functions and 
response functions can be defined as before, see eq. ( |5.2S| ). However, and in contrast to the non- 
equilibrium case, time translation invariance is expected to hold and we should have 

C = C(t-s;r) , R = R(t-s;r) (5.65) 

The response functions should transform covariantly under local scale transformations, as originally 
proposed in ]Tj|. In consequence, carrying over the treatment of section 4, R is again given by 
eq. ( |5^2| ), viz. 

R( T ; r) = (0(r, r) J(0; 0)) = T - 2xi/z $ (rr- 1/z ) (5.66) 

but with the additional constraint x\ = X2 coming from time translation invariance (alternatively, 
this constraint may be written as Xr/z = 1 + a). The scaling function $(w) is given as before by 



eq. ( 4.16p . The constraints discussed above in section 5.2 for the non-equilibrium response also apply. 
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Since at equilibrium the fluctuation-dissipation theorem 

T c R(r;r) = -^-C(r;r) 



(5.67) 



holds [19, 71|, where r = t — s, the scaling form of the two-time correlator may be given as well. 
Integrating, we have 



C{t;r)=C{0;r)-T c zr 



z-2x 



(5.68) 



Here C(0;r) = (<f)(t,r)<j)(t,0)) is the equal-time correlator, which is t- independent, because time 
translation invariance holds at the equilibrium critical point. 

We point out that the exponent x% which enters Q5.6S ) is not a static critical exponent. Rather, 
it is related to the static scaling dimension x± by 



2x{ 

,(g) 



(g) 



2x\ — z 



(5.69) 



For example, for a 2D critical point, the xf are the scaling dimensions which can be obtained from 
the representations of the Virasoro algebra ||. In terms of these, we have, with u = rt~ x l z 

R(t;r) = r 1 - 29 * /*^ (rt" 1 /*) (5.70) 

C(t;r) = C(0;r)-T c zr~ 2x ^ Tdww^-^iw) (5.71) 

J u 

Because of the known behaviour of $>(w) at the w — > and w — > oo boundaries, it is easy to see 
that the integral in ( |5.71| ) is convergent if x^ > and z > 0. Since we are at criticality, we expect 
C(0;r) 



,-2x 



1 and finally obtain 



C(t; r) = r- 2 ^ (c + C x £ dw w 2x ^- l ^{w)j , u = rf 



l/z 



(5.72) 



where Co,i are non-universal normalization constants. 

We now compare the scaling form ( |5.70|) with the prediction of dynamical conformal invariance for 
a non-conserved order parameter R(t,r) ~ t _1_2:c ^ exp(— r^/t) , up to suppressed non-universal 
constants. In our theory, we have found a simple exponential scaling function in two cases: (i) z — 2 
with /3i+7i = A^/2^0 and (ii) z — 1 with (3\ + 71 = 0. In these special cases, the two-time 
correlator takes the following form: in the case z = 2, M. ^ 



C7(f,r) 

and in the case z = 1, /3i + 71 = 



-2x- 



,(g) 



Co + C l7 (4 gJ ,— y 



C(t, r) = r" 2 *i s) [C + d 7 ( 2xS s) , 2 7l ^ 



(5.73) 



(5.74) 



respectively. Here y(a,z) is an incomplete gamma function, see eqs. (6.5.4), (6.5.29) in [0] and Co,i 
are normalization constants. In all other cases, the scaling functions will take a form quite distinct 
from these two examples, see figures [l] and [| We recall that the fundamental hypothesis of dynamical 
conformal invariance was conformal invariance in space []TS| and not in time, as we have assumed 
throughout this paper. 

Tests of the consequences (|5.70Jk72 ) of the hypothesis of local scale invariance in equilibrium 
critical dynamics in specific models would be most welcome. 
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6 Conclusions 



Our starting point has been the heuristic idea that it might be possible to extend strongly anisotropic 
or dynamical scaling from mere dilatation invariance with a given anisotropy exponent 9 (or dy- 
namical exponent z) to a larger dynamical symmetry involving local scale transformations with a 
space-time-dependent dilatation factor. We have studied this idea by seeking confirmation on the 
purely phenomenological level of being able to reproduce certain two-point function in the context 
of specific models. The agreement found provides evidence in favour of, but does not prove, this 
hypothesis of local scale invariance. 

In attempting to construct local scale transformations for an arbitrary anisotropy exponent 9, 
we have tried to follow the two known cases of conformal and Schrodinger invariance as closely as 
possible. The common feature of these groups is the presence of conformal transformations in time 
and we have made this the central assumption in our formulation of local scale invariance. 

In carrying out the construction of infinitesimal local scale transformations, we have seen that 
there exist two types, the first one (Typ I) being related to strongly anisotropic equilibrium systems, 
while the second one (Typ II) describes systems with time- dependent dynamical scaling. The main 
properties of local scale invariance when applied to quasiprimary operators are collected in table |2|. 
Conformal invariance and Schrodinger invariance are recovered as special cases, for 9 — 1 and 6 = 2, 
respectively. On the other hand, if 9 ^ 1,2, the generators only close into a Lie algebra on certain 
states only and thus form a weak Lie algebra. 

Local scale transformations form a (weak) dynamic symmetry group of the equation of motion 
Sip = of certain, in general non-local, free-field theories (see table |2] for the precise form of S). 

In addition, two-point functions formed from quasiprimary operators satisfy certain linear (frac- 
tional) differential equations, from which the form of these two-point functions can be determined. 

Indeed, these explicit predictions (see eqs. ( |4.23| , [i~2"5"|) for Typ I with N = 2/9 6 N and eq. (|4.16|) 



for Typ II) allow for a test of the applicability of our notion of local scale invariance in concrete 
models. These tests have been performed in two different settings: 

1. Uniaxial Lifshitz points are a classic example of a strongly anisotropic equilibrium system. The 
form of the spin-spin correlator at the first-order Lifshitz points in both the ANNNS model 
and the 3D ANNNI model and also at the second-order Lifshitz point in the ANNNS model 
agrees fully with local scale invariance, for 9 ~ | and 9 = |, respectively |55|, |90 . 



2. Dynamical scaling occurs in the aging behaviour of simple ferromagnets which undegoes phase- 
ordering kinetics or non-equilibrium critical dynamics. We found a particularly simple scaling 
form of the two-time autoresponse function in the aging regime 



= s-'-°f R (t/ S ) 

h=0 



5h(s) 

f R (x) = r x A (x-l) B (6.1) 

where the values of the exponents A and B can be matched to known dynamical exponents, 
see section 5.2. This form can be reproduced in a large variety of models, with the dynamical 
exponent z taking values both below and above 2, notably in the 2D and 3D Glauber-Ising 
models, the 0(n)-symmetric kinetic model A to first order in the e-expansion, several variants 



of the kinetic spherical model with non-conserved order parameter and brownian motion |58 
|I§ 0, 
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The form ( |6.1|) for the scaling function fn(x) is likely to be very robust and of broad validity, 
see section 5.2. 



We stress that the 3D ANNNI model and the 2D and 3D Glauber-Ising models cannot be expressed as 
free-field theories. This renders the confirmation of local scale invariance in these models non-trivial. 

Technically, our results depend on the construction of commuting fractional derivatives. That 
property was essential in the derivation of the fractionall differential equations satisfied by the scaling 
function of the two-point function. However, these differential equations themselves do not depend 
on any other property of the fractional derivatives. On the other hand, a specific choice must be 
made for the fractional derivative if an explicit solution is requested. At present, there is not yet any 
test available which would inform us whether the specific form used in this work or else any of the 
known alternatives is realized in the context of local scale invariance. 

The possibility of local scale invariance has recently been discussed (under the name of repara- 
metrization invariance) in the aging of certain glassy systems, where the covariance of correlators 
and response functions was studied [^, ^]. In the models studied there, the scaling dimensions 



= of the order parameter and it = 1 of the response operator take rather simple values, 
however. The covariance of non-equilibrium two-time correlators under time reparametrizations also 
plays a role in a recent study on general constraints on the scaling of these [173] . Finally, a recently 
introduced model for a fluctuating interface with z = 1 can be described in the thermodynamic limit 
in terms of the characters of a conformal field theory 



All in all, we have proposed a new type of dynamical symmetry which might generalize usual 
anisotropic or dynamical scaling. The existing phenomenological confirmations provide evidence for 
local scale invariance to be realized, at least to a very good approximation, in certain statistical 
systems. Of course, the theory must be built further and in particular, one would like to be able to 
show that the field theories underlying these scale-invariant statistical systems indeed satisfy local 
scale invariance. Work along several of these lines is in progress. 



Note added in proof: 

For the critical weakly disordered kinetic Ising model with a non-conserved order parameter, Calabre- 
se and Gambassi ||118|| find, to one-loop order, the response function at vanishing external momentum 



q to be in full agreement with the prediction of local scale invariance. 
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Appendix A. On fractional derivatives 



Fractional derivation and integration is a well-established topic, see [99, 79, 91, 60] for reviews. 
However, for a derivative operator d a of real order a the commutativity d a d b = d a+b = d b d a is 
not trivial. The well-known Riemann-Liouville and Griinwald-Letnikov fractional derivatives do not 
satisfy it. Since our construction of local scaling operators is built around this property, we shall 
present here fractional derivatives in a self-contained manner such that commutativity is guaranteed, 
along with some other simple calculational rules. In particular, the definition used here allows for 
the solving of fractional differential equations of rational order by series expansion methods. For 
simplicity, we merely consider functions f(r) of a single variable r. 

Consider a set E of numbers e such that any e G E is not a negative integer, e ^ -(n + 1) with 
nGN. We call such a set an E-set. Let I be some (possibly infinite) positive real interval. We define 
the Ai -space of generalized functions associated with the E-set E 



M:=M E (I,R) = {f:IC 



f(r) = J2 fer e + Yl F ^ [n) (r)-J e eR,E n eR\ (Al) 



e£E n=0 



where 5^ n \r) is the n th derivative of the Dirac delta function. The part of / G M. parametrized 
by the constants f e is called the regular part of / and the part parametrized by the E n is call the 
singular part of /. A well-known theorem p. 81] states that any generalized function /(r) 
concentrated at r = is indeed given by a finite sum ^2 n F n 8^(r). It is understood throughout 
that only finitely many of the coefficients F n in (|Al]) are non-vanishing. Furthermore, we assume 
that the regular part in ( |A1| ) converges 'well enough' that all 'reasonable' operations can be carried 
out (see Lemma 2 below for a sufficient condition). For example, we could choose E — N and take 
Ai to be the set of analytic functions on I. More generally, if we take E = fiN + A with \i > and 
A 7^ —{^{n + 1) + m + 1) with n,m e N, Mb is the space of functions of the form r x f (r M ) with 
f(r) analytic. Although everything here is specified in terms of real numbers, the formal extension 
to complex-valued functions is immediate. 

Definition: Let a G M, E be an E-set and let E' := {e'|e' = e — a; e G E}. Analogously to (|Alj ) one 
has the space M.' = M. E i. An operator d a : M. — > M.' is called a derivative of order a, iff it satisfies 
the properties: 

i) d a (A/(r) + ng(r)) = Xd a f(r) + /jd a g(r) VA,/iGl and all f,geM 

") d\ e = r , r(e + ! ) 1 / e - a + 5 a , e+n+1 T(e + l)tfW(r) (A2) 
r(e-o+l) ^ 

„— 1— n— a 00 

iii) d a 5 W (r) = _ + ]T <5 a , m 5(" +m ) (r ) 

v y m=0 

where T(x) is the Gamma function. In particular, it follows from (A~2) that the prefactor for any 



monomial r n 1 with nGN indeed vanishes. For our purposes, we may consider the set E' therefore 
also as an E-set and Ai' as an .M-space. 

In particular, it is implied that d a can be applied term-by-term to any function / G M. Often, 
we shall also write d a = 9" if we want to specify explicitly the variable r on which d a is supposed to 
act. We point out that d a is not defined on negative integer powers r _n ~ 1 with n > 0. 

In order to show that this definition is not empty, recall the definition of fractional derivatives as 
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given by Gelfand and Shilov 



PI 



p. 115]. For the real line, they consider the generalized function 



r > 
r < 



(A3) 



and for generalized functions concentrated on the half-line r > 0, they define 

-a-l 



d a f(r) := f(r) * 



1 



F -a) r -a) 



dp f(p) (r - p) 



-a-l 



(A4) 



o 



It is understood that the integral must be regularized ||42|| . The conditions (|A2| ) are immediately 
verified [42], using J^dtt a ^ 1 (l — t) b ^ 1 = T(a)T(b)/T(a + b) and analytic continuation. For comparison 



with the Riemann-Liouville/Griinwald-Letnikov/Marchaud fractional derivatives D a = oD®, we recall 
that D a r e = (r(e+l)/T(e— a+1)) r e_a p9|, ff9|, |9l|, |60| . The singular terms which arise in the definition 
are absent. 



The practical interest of the definition ( [A2] ) comes from the following simple rules for calculation. 

Lemma 1: If E is an E-set, M. the associated Ai-space, f G Ai and d a the derivative of order a 
such that all d a f considered below exist, one has on At 



d a+b f(r) 
[d a ,r]f(r) 
d?f(\r) 



d a d b f(r) = d b d a f(r) 
ad^fir) 
X a d a Xr f(Xr) 
\- a d a r f{r) 



where A > is a real constant. If f is analytic without singular terms and g G Ai, one has 



d a (f(r)g(r)) = J2 



d e f(r) 
dr e 



d a -"g{r) 



(A5) 
(A6) 
(A7) 
(A8) 



(A9) 



where d e /dr e are ordinary derivatives of integer order I. These rules are the natural generalizations of 
the familiar properties of the usual derivative. The commutator [., .] is defined as usual. Well-known 
fractional derivatives such as the Riemann-Liouville, Grundwald-Letnikov or Marchaud fractional 
derivatives satisfy the commutativity relation ( |A5| ) only if further conditions are imposed on a and 
b or on the function /(r) f99, TS, 31}| (see also the example (|A14|) below). On the other hand, (|A5|) 
does hold for the Gelfand-Shilov definition ( |A4[ ) |4"2] , pi| ] , the Weyl fractional derivative |7^] and for 
the recent complex multivalued definition presented in [|115|| . The generalized Leibniz rule ( |A9| ) is 
usually proven for both /, g analytic or infinitely often differentiable |^T| . In J79| , / is analytic 
and g G M. m, with the E-set M = N + p, p > — 1 (logarithms are also admitted). Generalized (and 
simpler) Leibniz rules based on the convolution product were discussed in |70| . 

So far, the definition of d a was performed for variables r > (and we should have written r e + 
instead of r e everywhere) . We may use ([AT]) with A negative to formally extend the definition of d a 
to any value of r ^ 0. Then 



AG,Aq) remain valid. 



Proof: We show that the rules of the lemma can be reduced to the basic properties 
to prove (|A5|) , we consider separately two cases. Let 



ind 



if b^N 
if ben 



2|). In order 



(A10) 
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First, consider 

9 a d b S(r) = d a (f^- + <W ( 

V ^ ' meN 



m) 



r 



-1— b— a \ I r^—l—m—a 



meN 



v ' \ v y neN / meN \ v ' neN 

and this also implies commutativity on 5^ n ^(r) = <9 n 5(r). Second, we have 
d a d b r e = d a ^ -I^±^ r e - 6 + g^, e+ » +1 r(e + l)^)(r)j 

U — 1,af r(e-6 + l) ^r(e-a-6+ 1) 

+Vw ( w (e + ^ " 1 ~ n " a + E ^ r ( e + i)^ (n+m) ( 

neN \ v 7 meN 

= r( e -t-a + i) r " t - + £ + U^M 

v y meN 

= d a+b r e 

Having thus checked the claim on the 'basis set' spanned by r e and 5^ n ' it holds on Ai by linear 
superposition. Clearly d^.d b r j = d h r d^j . In the sequel we shall need the identities 

r5 {n \r) = -n<5 (n - 1} (r) , 8 ( ~ n \\r) = A _ri ~ 1 5 (n) (r) (All) 

To prove (|A6|) , consider 

[d r V] r e = d?r e+1 - rd?r e 



T(e + 2) r( e + i) 



r e-a+l 



T(e-a + 2) r(e-a + l) 

oo 

+ E 5 a , e+n+2 r(e + 2)5( n )(r) - E <W+ir(e + l)n5< n >(r 

n=0 n=0 



oo 



n=0 
aa— l„e 



r(e+l-(a-l)) 
= ad^V 

where in the third line the first identity ([All ) was used. Next, 
[<9 a ,r]5 {p) (r) = -pd a 5 i - p ~ 1 \r) -rd a 5 (p) {r) 

\T(-a-p + l) Y{-a-p)J 



_^ r -l-p-(a-l) 



m=0 

oo 



r(-(a-i)-p) 
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where the first identity ( |A11| ) was used again. Having checked ( |X6|) on r e and S^ n \ it holds on M. 
by linear superposition. For (|A7|), we set s = Ar and use the second identity ( |A11| ). Then 

/ oo \ / oo 



n=0 



n=0 



J]A e / e r(e+l) 

oo 

+ ^ J F 1 „A— 1 



r 



r 



r(e — a + 

1— n— a 



m=0 



n=0 



T(—a — n 



+ Ys<rn5 [n+m) {l 



m=0 



= AM ^/ e 9t(Ar) e + ^F n 9 A V (n) (Ar) =X a d:f(s) 

\e£E n=0 / 

For the forth relation flA"8|), let again s = Ar and, using 

^r/(r) = d^fiX-h) = A _a 9"-i s /(A _1 s) = A- a 9 r a /(r) 



The fifth rule (|A9|) may be obtained in the same manner, but we shall prove below it with the help 
of the commutator identity ( |A17| ) which permits a shorter proof. q.e.d. 

In particular, one can set g(r) = 1 in eq. ( |A9| ) . Then, for / analytic and non-singular one has the 
regular part, see also |)9], [79], (Jl], ITT 



d a r f{r) 



ireg 



E 

£=0 



d e f(r) 



T(£-a + l) dr £ 



-E 



r(a + 1) sin(7r(a - £)) e d e f(r) 



(A12) 

which expresses the regular part of <9 a / in terms of ordinary derivatives. Clearly, if a — > k E N, one 
recovers from the second form ( |A12| ) via d^f(r)\ — > d fc /(r)/dr fc the ordinary derivative of integer 



order > 0. On the other hand, if a = k + a with k EN and < a < 1, we have 

^n—k—a 



<t +Q exp(r) = ^ 



r- fc - a ^a,i_ fc _ a (r) 



r(n+ 1 — — a) ^ r(n + 1 - a) 

ra=0 v ' n=—k v ' 



(A13) 



where E a ^(r) = Y^k=o z /r(atA; + /3) is a Mittag-Lefner type function. If a = 0, we recover 
c^f exp(r) | reg = exp(r), but that property does not hold anymore if a > (for the Weyl fractional 
derivative one has d a w^ r = e r f7"9fl). 

Finally, it may be useful to illustrate the commutativity of <9 a on some example. Following Miller 
and Ross [|7^, p. 210], take a positive integer q and let 



q-l 



(A14) 



k=0 



For the Riemann-Liouville derivative D a = qD^, one has indeed D 1 / q z(t) = at(t), but D 2 / q and 
are clearly different, since 

D 2 ^e(t) = a\(t) + ^ « 2 e(t) = D 1 ^ (L>^ e (0) 
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and this fact had led Miller and Ross to the definition of sequential fractional derivatives [79 j. On 
the other hand, using ( |A2| ) , we find 

d] /q t{t) = at(t) + 5(t) 

d 2 t 'h{t) = ah{t) + + a8(t) = (^e(t)) 

and the singular terms are seen to be essential for the commutativity property (|A5|). 

Consider an E-set E which is countable and ordered. Therefore, the elements e G E can be 
labelled by an n G N, viz. e = e n . Let v n := e n+ i — e n > 0. Call such an E-set E well- separated 
with separation constant e, if there is an e > such that v n > e. For the regular part of a function 
/ G M.E-, we have 

oo 

/( r )lrcg — ^ j fe r ~ ^ ] fn r i fn ■ /e n 

eeB n=0 

Questions of existence of d%f{r) and its relation to ordinary derivatives are dealt with in the following 

Lemma 2: Let E be a well- separated E-set with separation constant e, f G M.e, e n > ; e n — a > 
and 

p- 1 := limsup \f n \ 1/e " > (A15) 

n— »oo 

T/ien i/ie following holds. 

(i) f{r) converges absolutely for \r\ < p. 

(ii) If u n /e n < B for some constant B, d£f{r) converges absolutely for |r| < pmin (l, (1 + B)~ a / £ ) . 
(Hi) If f : I — > R is analytic with a radius of convergence p > around r = 0, then the series 
( \A1Q) for 9"/(r) converges absolutely for \r\ < p/2. Property (i) is well-known for psi-series and 
still holds if only e n / Inn — > oo as n — > oo ||62|| . The conditions imposed here are sufficiently wide 
to include functions of the form r x f (r^) with /(r) analytic, /i > and A ^ —pan — n, n, m G N, 
which is enough for the applications we have in mind. In this case, we have effectively B = 0, since 
v n j e n = pi/ {pin + A) — ► as n — > oo. 

Proof: The conditions e n > and e n — a > are only needed in order to make the regular parts 
of /(0) and d a f(0) well-defined. Since the singular parts of / and d a f are finite sums, they do not 
affect the convergence and can be suppressed here. 

(i) In order to show the convergence of f(r), consider iV G N sufficiently large. We then have the 
remainder stimate 

oo oo oo 

Rn := £ |Ar*| = £ l/V-rp < ^ ^ 

n=N n=N n=N " 

which holds for iV sufficiently large.0 Since is well-separated, we have e n — > (n — A^)e and 
e^v > e + Ne. Therefore, for |r| < p 



Rn < 



E 

n=7V 



-ejv 



< 



E 

n=7V 



e(n-AT) 



< 



eN 



\r/p\ e 



1 - |r/p| £ 



and -Rat — > as N — > oo. 

(ii) We need the asymptotic identity (6.1.47) in for z — > oo and constants a, 6 



T(z + a) 
T(z + 6) 



a—b 



{l + Oiz- 1 )) 



(A16) 



13 



More precisely, let 5 > and consider \r\ < p — 5. 5 can be made arbitrarily small if N is large enough. 
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For the convergence of 



<9 a /(r)| rcg = X> 



n=0 



T(e n + 1) 
r(e n - a + 1)' 



we consider the remainder 

oo 

Qn:=J2 



n=N 



/, „ r(en + 1 L r e -° 



r(e„ - a + 1) 



-IE 

n=JV 



r(e n + r 



r(e n - a + 1; 



|/Ve„ r |-<| r -a|^ e a 



where eq. ( |A16| ) was used and N was taken to be sufficiently large. Now, we use the following known 
fact P5| : if y n > and X^Lo^« < 00 an d furthermore \x n+ i/x n \ < y n+ i/y n , then the series J2^Lo x n 
is absolutely convergent. We apply this to the sequence x n := e^\r/p\ £n and have the estimate 



x n+l 


(r) 


r 




(" 




a 

) 


r 






P 








P 



< {1 + BY 



'n+l — e-n 



Therefore, if we take y n — (1 + B) an \r / p\ £n , it is only left to prove that Yl^LoVn is convergent. But 
this is obvious, since for |r| < p, one has 



oo oo 



Q' N :=Y,yn = Y,^ + B y 



n=N n=N 



< 



£ (1 + B) ° 

n=iV v r 7 



which indeed tends to zero for iV — > oo, if (1 + B) a {\r\/ p) e < 1. 

(iii) The series f(2r) = X]fcLo(^')~ lrfe f <yk \ r )i where f^ k \r) := d fc /(r)/dr fc , converges absolutely for 
|2r| < p, thus 

r j("+l)( r ) 



< 1 ; if \r\ < p/2 



n + 1 /(»)(r) 

If a e N, there is nothing to show. If a ^ N, let b = —a and recall that for a G 



a\ a(a — 1) . . . (a — £ + 1) 
= I\ 



,b(b + l)...(b + £-l) 
£\ 



We have from the first form of ( |A12| ) the estimate 



<^/(r)L g = d; b f(r) 



Ircg 



< r° 



E 



b(b+l)...(b + £-l ; j m 



£\ r(£ + b) 



oo 
1=0 



We call the £* coefficient in this series pi and have (since b is not a negative integer) 



P-n+l 



(b + n)r ( r ) 



< 1 ; if \r\ < p/2 



p n (n + l)(n + b) /(»)(r) 

which implies absolute convergence. 

Practical calculations with d a are simplified by the following 
Lemma 3: (i) If d a is the derivative of real order a and n e N, one has 



q.e.d. 



k=i ^ ' 



(A17) 
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(ii) If f is analytic without any singular terms and g G M. the relations ( \Ab}\A£\ , \An[ ) are equivalent. 
(Hi) If f G M. but without singular terms and n G N, one has with the ordinary derivative p*\x) = 
d e f(x)/dx e of integer order £ 



[f(dr),r n )=J2{ n £ ) rn ~ e f {i) ( d r) 

Proof: (i) We proceed by induction over n. The case n — 1 is the identity 
step, let g e M and consider, using (|A6|) again twice 



(A18) 

For the induction 



[a,", g(t) = [3,°. «"] tg(t) + 1" [8? , t] g(t) 



E 

v fc=i 

' n 

\fc=l 

n 

+E 

fc=l 



jfe + 1 



77 



fc! t n+1 ~ fc ar fe 



(fc + 1)! r +1 -( fe+1 ) ar (fc+1) + arar 1 



Ln + lror 1 + it( a k ) 

\ k=2 ^ ' 

n+1 

E 



n + 1 
k 



kit 



/c aa—k 



ga-k + 



n + 1 



{n + l)\dr n - l \g{t) 



k=l 



n + 1 
k 



k\ t n+l ~ k d?~ k g(t) 



and the assertion follows. 

(ii) For an analytic functions without singular terms, one has f(r) = Y^n=o fn rn an d f° r all A; G N 
the ordinary derivative of t n is 

d k f n 



dt k 



kit 



n—k 



From (|A17 ) we have 



9?(t n g(t)) = £ ( t ) ( I ) Mr- k ar"g(t) 



and (|A9|) indeed follows, under the stated assumptions on /. Conversely, starting from the generalized 
Leibniz rule (|A9[) and setting f{t) = t n with n G N we have 



[81 T] g(t) = % (t n g(t)) - T %g(t) = 



k=l 



k+n 



and we recover indeed ( |A17| ). The special case n = 1 then reproduces fl/Of). 
(iii) For / G .M non-singular one has /(r) = Y^ e / e r e . Therefore 



lf(dr), 



where in the second step ( jXq ) was used. The assertion now follows immediately by induction over 
n. q.e.d. 
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Appendix B. Generators of the Schrodinger algebra for d > 1 



We list the generators of the infinite-dimensional Lie algebra of the Schrodinger group in d = 2 
spatial dimensions. They read 

X n = -r +1 ft - ^nndr + r 2 d 2 ) - ^p^Mt^rl + r\) - |(n + l)t n 



y« = -t m+l ' 2 d 1 - (m + TWt" 1 - 1 



/2 n 



M n = -t n Al 

where o^- = <9/ <9r,- with j = 1,2 and where .M is the mass. Here rj G R are the two spatial coordinates. 
The indices n G Z and m G Z + |. With respect to the d = 1 case treated in the text, there are 
now two sets of generators for generalized Galilei transformations and the new generator R of spatial 
rotations. 

A straightforward calculation gives the commutators = 1,2): 

[X n ,X n i] = {n — m)X n+n i 

[X n , M n >] = -n'M n+nl 

\?£\Y$\ = 6 hJ (m-m')M m+m , (B2) 

[Y%\M n ] = [M n ,M„,] = 

[X n ,R] = [M ni R] = 

\X m \ R] yJ^ ) 
[^m ^ ) -R] yJ^ ^ 

which closes into an infinite-dimensional Lie algebra. Similarly, the extension coming from including 
a parameter B 2 o in the generators (see ( |3.28| , [3T34f) for the ID case) can be written down straightfor- 
wardly. 

The special case of the finite-dimensional Lie subalgebra, corresponding to ( [2.15| ) for d — 2, is 
given by the set 

{X_i A l, yli/2,+l/2> ^-1/2,4-1/2' Mo, R}- (B3) 

The generalization to spatial dimensions d > 2 proceeds along the same lines. 



Appendix C. Infinitesimal generators for generalized Galilei 
transformations 

As an alternative to the construction of local scale transformation given in the text, we present here 
a different construction of the infinitesimal local scaling generators and proceed to the calculation of 
scaling functions. This had been the first case where generators of a local scale invariance could be 
explicitly constructed and scaling functions could be found. 
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Starting from the generators 



X- X = -d t , X = -tdt - ~rd r , K_ 1/2 = -d r (CI) 
which satisfy the commutation relations 

[X ,X- 1 ]=X- 1 , [X_!,y_ 1/2 ] =0 , [X ,y_i /a ] = ^-1/2 (C2) 

we want to find, for 9 as general as possible, a generalized Galilei transformation Yx/ 2 such that a 
closed Lie algebra results. In particular, we want to construct this algebra A such that Xq acts as 
a counting operator, that is for any generator A e A, we have [Xq, A] = a A. Motivated from the 
Schrodinger case 9 = 2, where Yx/2 = —td r — Air (and after having tried out many different forms), 
we make the ansatz 

Yx/2 = -td r -M(d r )r (C3) 
In what follows, we use the formal properties of the derivative d°: as defined in appendix A. 
First, we formally calculate the commutator, using ( |A18| ) 



[X ,Y 1/2 ] = 6 —^td r - l -M\d r )d r r + ~M{d r )r = - d —^Yx/2 (C4) 

where the last equation is motivated from the first term in the generator Yx/2 an d our construction 
principle. This leads to 

(9 - l)M(d r )r = -M'(d r )d r r + M(d r )r (C5) 
If we let x := d r , we find (9 — 2)M.(x) = —Ai'(x)x which has the solution 

M{x) = M x 2 - 6 (C6) 

where Aio is a constant. Next, we find the commutators 

[x_ x , y 1/2 ] = -y_i /2 , [Yx/2, y-i/ 2 ] = -M d 2 r ~ e =■. m q (C7) 

where M = — Aiod 2 ~ e is a new generator. For 9 = 2, we recover the Galilei algebra, where M is 
central. Otherwise, we find through a formal calculation 

[X-i, M ] = [y_i /2 , M ] = , [X , M ] = -^M , [Yx/2, M ] = -(2 - 9)M 2 d^ 28 (C8) 

and we see that we must define a new generator iV. In the special case 9 = 3/2, we have iV := —^A4q. 
Then 

[y 1/2 ,M ]=iV (C9) 
and is central. Therefore, if 9 = 3/2, the set 

A := X , Y_x/2, Yx/2, M , iV} (CIO) 

closes as a Lie algebra and satisfies the condition that Xq acts as a counting operator. 

We now derive the form of the two-point function covariant under the transformations generated 
by the set A. Scaling operators will be characterized by their scaling dimension x and their 'mass'. 
For the Schrodinger invariant case 9 = 2, scaling operators are doubletts {(f), <$*), with 'masses' 
(A4, — M.), where M. is a non-negative constant |]53||. Here, for 9 = 3/2, it turns out that scaling 
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operators are quadrupletts 4>^ a \ with a = 0,1,2,3 and 'masses' Wl := i a M, where again M is a 
positive constant. We consider two-point functions of the form 



F = F^\h, t 2 - ri) r 2 ) = (0i Q) (ti, r x )0f ] (t 2 , r 2 )) 

and the covariance conditions are 

2x 

X F = —F , X^F = Y± 1/2 F = MqF = NF = 
3 

where x = X\ + x 2 . Here, and K_i/2 are given in (|CT|), while the other generators read 

Y 1/2 = -td r - Mdl /2 r , M = -dRd l J 2 , N = -]m 2 



(Cll) 



(C12) 



(C13) 



Spatio-temporal translation invariance yields F = F(t,r), where t — t\ — t 2 and r = r\ — r 2 . 
Invariance under the action of N leads to the condition M\ = (— l)P~ a+1 Ml or alternatively 



iP - a +i Mi ^ where ^ = a + 1 mod 2 



(C14) 



Then the action of Mq on F is, using translation invariance and ( |A7|) 

M F(t,r) = -(M 1 i a d¥ 2 + M 2 i p+l dy 2 )F(t,r) 
= -i a (Mi + M 2 i p - a+l ) dl /2 F(t, r) 
= -i a Mt (1 - 9 r 1/2 F(t, r) = 



(C15) 



because of ( |U14 ) and thus F is always invariant under the action of M . To calculate the action of 
of Y1/2, we recall from the Leibniz rule ( |X9| ) that for f £ Me 



dl' 2 (r/(r)) = rdy 2 /(r) + iflrVa/( r ) 



(C16) 



and we find, using again (|A7|) 



- m x ( + i^/ 2 J - m 2 Udli 2 + \d^ 2 



F(t,r) 



-td r - Mii a (n - r 2 i 2 ^- Q+1 )) d)' 2 - -M{? (1 + i 2 ^-^ 1 )) d; 1 ' 2 



F(t,r) 



= -(td r + Tl(rdl /2 + d; 1/2 ))F(t, 
where we have set 9Jt := M\i a . Finally, 



X F(t,r) 



td t + -rd r ) F{t, r) 



(C17) 



(C18) 



Now, the scaling ansatz 

F(t,r) =r 2x/3 ^(u) , u = rr 2/3 (C19) 

solves the first of the covariance conditions flC12|) , while the last remaining condition ( C12|) leads, 
via (pT7) and (|A7|), to a fractional differential equation for the scaling function ty(u) 

(du + mud^ 2 + md~ 1/2 )^(u) = o (C20) 
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which coincides with eq. (|3.95| ) found for Typ Ha with /3 — and 7 = |07t. Inspection of the terms 
present in this equation leads to the following series ansatz 



(C21) 



n=0 

which promptly gives s = and the recursion relation, valid for all n > 1 

,'3n — 1\ „ /3n\ T 



(C22) 



This is best solved by rewriting it in the form ^ n+2 = ^ n (207t 2 /3)(n + 1) ((n + 2/3) (n + 4/3)) 1 and 
we finally have for desired scaling function 



071 



3/2 



3 ^ r (l + I) r (? + I) V 

, F2 U, I i, H, ™V) - ,/l^7 2 F 2 C 1, ^ 
" V 2' 3' 3' 3 / V vr V 6 6' 3 



(C23) 



where 2 F 2 is a generalized hypergeometric function and \l/o = ^(0) is the initial value. 

For a physical interpretation, recall that the four scaling operators (p^ had the 'masses' 071 = 
i a Mi, where Mi > and a = 0, 1,2,3. In addition, considering global scale invariance with either 
t = or r = 0, it follows that for u — > 00, the boundary condition ty(u) — > must be satisfied. 
That is possible in two cases: (i) when 071 > and (ii) when 07t is imaginary and the real part of the 
two-point function is retained. In these cases, the two-point functions read 



F (0,1) (t, r) = r 2x/3 ^(n) ; where 071 = Mi = M 2 > 
F, = l(FlV\t,r)+FM\t,r))=r*''^ 

; where OJt 2 = -Ml Mi = M 2 > 
The scaling functions ^i l2 (u) obtained form these are shown in figure BL 



(C24) 



Appendix D. 

We discuss the solution of the fractional-order differential equation ( |4.22| ) in terms of series expansion, 



following standard lines [|79| |91fl . As in section 4, we have N = Nq+ p/q with N E N. We make the 
ansatz 



n=0 



(Dl) 



and from substitution into eq. (|4.22j ), we find for N > 2 
aif (~{n + p + qN ) + s + 1 



n=—p—qNo 



n/<y _ 



V f -(n + pxi + gA^ xi) + s j a n W 9 = 



(D2) 
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-0.5 



Figure 9: Scaling functions ^1,2(1*) for the two-point functions Fj(t,r) = t 2x ^ 3 ^i(u), for M. = 1 and 
*o = 1- 



where the definition (|A2|) has been used. The singular terms can be dropped if v is positive. Since 
this equation holds for all v > 0, we obtain the following conditions. First, since ao 7^ 0, we must 
have 

a-iT (s + 1) .„ , 

7 — - — v = ( D3 ) 
r(-|-iVo + 2 + s J 

Second, we have = for £ — 1, . . . ,p + qNo — 1. Finally, we get the recurrence 

ai T U(n + p + qN ) + s + l) /j x 

7 r y -a n+p+qNo = [~(n + pxi + qNoXx) + s a n (D4) 

r f 1 + 2 + sj v<? / 

From the first condition eq. ( |L)3| ) , we find the possible values of s, namely 

s = s m := - + m , m = 0, 1, . . . , N — 2 (D5) 

q 

Negative values of s would lead to a singularity as v — > and are therefore excluded. The second 
and third conditions are solved by writing n = {p + qN )£ and by letting ut := a^^^t = a n- We 
then find 

n 2 I, , s m + c \ /. , \ + C\ r(w + »„ + i) 
" ( f + —J + —J r W < + 1) + s m + 1) (D6) 

The AT — 1 linearly independent solutions are, with e = N — N Q 

00 

«(m)(«) = J>^ +£+m 



£=0 

£+mu r(z + ( Sm + Q/iv)r(/ + (s m + 1)/N) T(s m + 1) 

7 r(( Sm + c)/iv)r(( Sm + i)/iv) r(iw + Sm + 1) V «i ' 



where m — 0, 1, . . . , Nq — 2 and cuo.m are arbitrary constants. It is easy to see that these series have an 
infinite radius of convergence, provided N > 2. Although this solution was only derived for rational 
values of N = N Q + e, we can make the analytical continuation to all real values of N. 
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We have seen in section 4 that there are solutions of eq. ( |4.22| ) such that Q(v ) ~ for v large. 
These may be constructed from the series ( |L)7| ) in the same way as done for N integer in section 4. 
On the other hand, for v small, one has Q^(v) ~ v £ . Therefore, the boundary condition Q(0) = 1 
cannot be satisfied for e > 0, that is any non-integer value of N. 

This difficulty with the v — > boundary condition is a specific property of the fractional deriva- 
tive ( |A2|) which it has in common with the Riemann-Liouville derivative. It is known that initial 
conditions for ordinary fractional differential equations are specified in terms of fractional derivatives 
d a f(0) [79, |99], [H], |6(J and not in terms of ordinary derivatives f^ N \0) of integer order N. Indeed, it 
has been suggested to avoid this problem by using the fractional Caputo derivative instead ]91 , |60 |. 
However, the Caputo derivative does not commute and it appears to be an open mathematical 
problem if the Caputo definition can be modified such as to obtain a commutative operator. 
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